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Abstract 



We study the no gravity limit Gjv — * of the Ponzano-Regge amplitudes with massive 
particles and show that we recover in this limit Feynman graph amplitudes (with Hadamard 
propagator) expressed as an abelian spin foam model. We show how the Gn expansion of the 
Ponzano-Regge amplitudes can be resummed. This leads to the conclusion that the effective 
dynamics of quantum particles coupled to quantum 3d gravity can be expressed in terms of 
an effective new non commutative field theory which respects the principles of doubly special 
relativity. We discuss the construction of Lorentzian spin foam models including Feynman 
■ propagators. 
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1 Introduction 



Spin Foam models offer a rigorous framework implementing a path integral for quantum gravity [1]. 
They provide a definition of a quantum spacetime in purely algebraic and combinatorial terms and 
describe it as generalized two-dimensional Feynman diagrams with degrees of freedom propagating 
along surfaces. Since these models were introduced, the most pressing issue has been to under- 
stand their semi-classical limit, in order to check whether we effectively recover general relativity 
and quantum field theory as low energy regimes and in order to make physical and experimental 
predictions carrying a quantum gravity signature. A necessary ingredient of such an analysis is the 
inclusion of matter and particles in a setting which has been primarily constructed for pure gravity. 
On one hand, matter degrees of freedom allow to write physically relevant diffeomorphism invariant 
observables, which are needed to fully build and interpret the theory. On the other hand, ultimately, 
we would like to derive an effective theory describing the propagation of matter within a quantum 
geometry and extract quantum gravity corrections to scattering amplitudes and cross-sections. 

In the present work, we study these issues in the case of a three-dimensional spacetime. Three- 
dimensional gravity is seen as a toy model in the search of a consistent full quantum general 
relativity. It is known to be an integrable system carrying a finite number of degrees of freedom. 
Nevertheless, despite its apparent simplicity, we face many of the same mathematical and inter- 
pretational problems than when addressing the quantization of four-dimensional gravity. In this 
context, the Ponzano-Regge model was actually the first model of quantum gravity ever written [2]. 
It gives an explicit prescription for the scalar product of Euclidean three-dimensional gravity as a 
state sum over discretized geometries. As the simplest non-trivial spin foam model, it is the natu- 
ral arena to investigate the particle couplings to quantum geometry and the semi-classical regime 
of the theory. Recently, this model has been studied in much detail. The paper [3] has tackled 
particle insertions and shown how they can be understood as a partial gauge fixing of the state 
sum model. As a result, explicit quantum amplitudes for massive and spinning particles coupled to 
gravity were constructed. The work in [4] established an explicit link between the Ponzano-Regge 
quantum gravity and the traditional Chern-Simons quantization, and identified a K-deformation of 
the Poincare group as the relevant symmetry group of spin foam amplitudes. In this article we 
follow the same line of thoughts focusing on the relation between spin foam amplitudes and usual 
Feynman graph evaluations. 

The first step is to show that the usual Feynman graph evaluations (with the Hadamard prop- 
agator) for quantum field theory on a three-dimensional flat Euclidean spacetime can be identified 
through a duality transform as spin foam amplitudes of an abelian model. This construction can 
be extended to the non-abelian case, in which case we describe Feynman diagrams for particles 
propagating in 3d quantum gravity. The natural mass unit for particles is the Planck defined solely 
from the Newton constant for gravitation mp = 1/k = 1/AnG. 

We define a gravity- less limit of the full quantum gravity theory as the G — > limit. And we 
show that the Ponzano-Regge amplitudes constructed in [3] reduce to the amplitudes of the abelian 
model with the standard Feynman evaluations: as required, the classical limit of matter insertion 
is the standard quantum field theory. 

Then, at G ^ 0, the spin foam amplitudes are to be interpreted as providing the Feynman 
graph evaluation of particles coupled to quantum gravity. We study the perturbative G expansion 
of the spin foam amplitudes. Remarkably, this expansion can be re-summed and expressed as the 
Feynman graphs of a non-commutative braided quantum field theory with deformation parameter 
G, which thus describes the effective theory for matter in quantum gravity. Any deformed Poincare 
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theory usually suffers from a huge ambiguity [5] coming from what should be identify as the physical 
energy and momenta since the introduction of the Planck scale allows non-linear redefinitions. This 
ambiguity can also be understood as an ambiguity in the identification of the non-commutative 
space-time. Our work show that the Ponzano-Regge model naturally defines a star product and 
a duality between space and momenta, therefore no ambiguity remains once we identify quantum 
gravity as being responsible for the effective deformation of the Poincare symmetry. 

This realizes explicitly, for the first time from first principles, the now popular idea that quantum 
gravity will eventually lead to an effective non-commutative field theory incorporating the principle 
of Doubly special relativity [6]. 

When we also take into account a non vanishing cosmological constant A, 3d Euclidean quantum 
gravity is described by the Turaev-Viro model whose coupling to matter was investigated by Barrett 
et al. [7, 8]. It is based on the quantum group su q (2). A and G now define natural length and 
mass units, bounding the physical quantities by above and below. We define different classical 
limits, which correspond to sending some of these units to zero or infinity. Interestingly, they all 
correspond to the same limit of the deformation parameter q 1. However they define different 
limits of the Feynman graph evaluations, depending on which units we use to describe the mass and 
length of physical objects. We define a hyperbolic spin foam models corresponding to a space-time 
at A < 0. We study the gravity- less limit G — > of the Turaev-Viro model and of the hyperbolic 
state sum model. 

Finally, we tackle the issue of causality of the particle amplitudes. Indeed, in the Euclidean 
framework, we have dealt up to now with the Hadamard propagator. It is non-oriented in time and 
non-causal. In the Lorentzian spin foam model, it is possible to discuss the different propagators 
-Wightman, Feynman and Hadamard- and write down the corresponding amplitudes in the non- 
abelian set-up. This means that we now have the explicit proposal for the amplitudes of the 
non-commutative quantum field theory describing the effective theory of 3d Lorentzian quantum 
gravity. 



2 Spin Foams as the topological dual of Feynman diagrams: the 
abelian case 

We consider the following Feynman diagram amplitudes associated to an oriented graph T for a 
massive relativistic scalar field in the 3d Euclidean spacetime: 

7 r = / Y[dx v JjG me (x t(e) -x s(e) ), (1) 

•* v e 

where v and e label the vertices and the edges of the graph T, m e the mass of the particles living 
on the edge e and G m the propagator of the theory. The x correspond to the positions of the 
particles/vertices of the graph (they are 3- vectors). Here we take: 

G m (x) = 55*1. (2) 

Let us insist on the fact that we are using the Hadamard function, and not the Feynman propagator. 
So we should call Feynman- Hadamard evaluation, but we keep it to Feynman evaluation whenever 
there can not be any confusion. We will later discuss the use and the implementation of the 
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Feynman propagator in the spinfoam framework, when dealing with the Lorentzian version of the 
models. In the present section, we show how we can recast the Feynman evaluation as a spinfoam 
amplitude through a duality transform. 

Let us start by switching to the Fourier transform. For this purpose we use the Kirillov formula: 



G m (x) = m [ d 2 ne im ^ H , (3) 
Js' 2 

then we change the space of integration from S 2 to M 3 using: 

L** = ik?L* 9m - m) - (4) 

I T = ! JT«» / 11^ ^tm ~ II e "' V(i '""~ i '""'- < 5 > 

v e e e 

Integrating over the x v , we finally get: 

/n^n<(r*-E*). «» 



This leads to: 



The 5(\p\ — m) fixes the norm of the momenta to be the mass 1 and the 6(J2Pe) at the vertices 
express momentum conservation. 

Let us now go to the topological dual and express this amplitude as a Feynman evaluation on 
the dual graph. A similar duality transformation on Feynman graphs was performed in [9]. 

More precisely, let us embed the graph T in a surface S. This amounts to defining faces / 
(as sequences of edges) on the graph, each edge belonging to two faces. Then one can define the 
topological dual, the dual vertices v being the faces, the dual edges e linking the v being associated 
to the original edges e and the dual faces corresponding to the initial vertices v. 

We would like to solve the constraint on the dual faces / (i.e at the vertices v) imposed by the 
S(J2e£df e f(c)Pe) (where ej{e) is a sign recording the orientation of the edge e relatively to the face 
/). From the algebraic topology point of view, p is a field valued on the edges e, so it is a 1-form. 
Its derivate dp is a 2-form ie a field valued on faces /. More precisely: 

{dp) j = (df,p) = e/(e)pg. 

e&dj 

Therefore we are imposing that dp = 0. If the (co)homology group H 1 ^) = Hi(S) is trivial, all 
closed 1-forms are exact. Then there exists a field u valued on vertices v such that p = du i.e 
Ps = Ut(g) ~ u s{e)-> whith t(e), s(e) being the starting and terminal vertices of eb. More generally H 1 
is generated by the cycles of S. So if S ~ S 2 , everything is straightforward. But if the genus of S 
is g > 1, then we also need variables uq attached to the cycles of S. These cycles C are defined, 
similarly to faces, as sequences of edges e. And now the solution to d X are given by: 

Ps = u m - u s( e) + e c(e)u c , (7) 

CBS 

1 S(\p\ — m)/2m is equivalent to the usual 8([p\ 2 — m 2 ). 
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where ec(e) = ± depending on the orientation of the dual edge e along the cycle C. 

A simple counting argument could reproduce the same output. Indeed we have E (number of 
edges) variables p s and V — 1 (V number of vertices) constraints. So we would need E — (V — 1) = 
F — x + 1 variables to parametrize the solutions to this linear system, with x the Euler characteristic 
of the surface S. As we use variables associated to the faces (i.e the dual vertices), let us allow for 
a closure relation on these face variables, so we would need F — x + 2 variables. If S is orientable 
then, this gives F + 2g variables: one variable per face and one variable per cycle. The same holds 
for non-orientable surfaces with F + g variables. 

Finally the Feynman graph evaluation reads: 

f ^i\p I ) 
ir = / Yi d ™v II d ™ c II 1^ — ~ with Ps = n *( g ) ~ Us ( s ) + ^2 e c(e) u c- (8) 

•^ R3 v C cycles e CSe 

The reader can find examples for S = S 2 and S = T 2 in appendix. 

More generally, let us consider a triangulated surface S and the graph T embedded in the 
triangulation (but not necessarily covering the whole triangulation) . We could more generally 
consider a generic cellular decomposition of S without loss of generality. We attach variables to 
every face/ triangle, or equivalently to all dual vertices v, and to every cycle of S. Then the Feynman 
evaluation simply reads: 

/r = / 3 11 du * 11 dU C 11 4^ 11 S ( U m ~ U *(e)), 

JK v£S C cycles of S eer e eeS\T 

(9) 

where the 5 functions on the edges in S \ T allows to 'collapse' the triangulation of S down to the 
graph T. 

We now embed the Feynman diagram in the three-dimensional spacetime and recast the Feyn- 
man evaluation in terms of the three-dimensional objects. Let us start with a three-dimensional 
triangulation A (or more generally a cellular decomposition), on which the graph T is drawn. We 
now would like to express the Feynman amplitude in terms of variables living on the structures 
(vertices, edges, faces or 3-cells) of A of its dual. We remind that the (2-skeleton of the) topological 
dual of A is called the spin foam. 

Let us choose a surface S C A in which T is faithfully embedded 2 and rewrite the formula (9). 
S is a collection of faces / (triangles more exactly if A is strictly speaking a triangulation) to which 
we associate vectors itf. Equivalently, one can consider the dual edges e* transverse to the faces / 
and the corresponding variables u e * . Then one can directly re-write the Feynman evaluation as: 

T f TT A~" TT A- TT 6 (l^W ~ g *( e ) + ^C3e S c( e )uc\ ~ mj y-r _ _ . 

J r = / 3 11 du f 11 du C 11 4^ 11 *K(e) " «.(e)), 

JK fes c cycles of s eer e ee5\r 

(10) 

where t(e) and s(e) are the two faces of S adjacent to the edge e £ S, and where e G C means 
that the link e dual to the edge e on the surface S belongs to the cycle C. We give an example for 
S = S 2 , A = S 3 in appendix. 

Let us now show how such an amplitude arises from spin foam models. Once the three- 
dimensional triangulation A is chosen, one can define the partition function of the (topological) 

2 We require that the components of S when the graph F is removed are all homeomorphic to a 2d disk. 
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spin foam model corresponding to the quantization of an abelian BF theory for the gauge group 
R 3 on the dual of A. One associate group variables -here vectors- u e * to every dual oriented edge 
e* (or equivalently to every face of A) and one requires that the holonomies around all dual faces 
/* are trivial: 

z A b) = / n^n^i. (ii) 

where we have introduced the notation 

= Yl e H e *R*- 

In order to introduce particles, it is more convenient to expand the 6 function in the partition 
function using (3). Then using the fact that a dual face /* is literally an edge e € A, the spin foam 
amplitude reads: 

za= [ u du e * [ n^n^r^- ( i2 ) 

e*gA* + egfl e€d 

Usually, a spin foam amplitude is presented as a product of weights depending solely on the l e 
variables after integration over the holonomies u e * . Nevertheless, this first order formalism, keeping 
both u e * and l e variables, is necessary to describe the insertion of particles. More precisely, inserting 
a particle on a given edge e with mass m e corresponds to the observable: 

oS"»<™.) - s -f - <*) 

One can obtain this expression as the abelian limit of the Ponzano-Regge spinfoam model discussed 
in the next section whose particle insertions have been described in [3]. Then inserting a whole 
Feynman diagram corresponds to putting particles on all edges of a given graph T. Inserting the 
corresponding observable into the partition function gives: 

Z f\r, {me }) = i n du e * n 5{ ^i~ me) n (1 4 ) 

J e* eer e e 

Now one would like to show that the previous spin foam amplitude (14) reproduces the Feynman 
evaluation (10). The two expressions already look very similar. But we would like to be more exact 
and identify a particular 3d triangulation A (or more precisely a class of triangulations) such that 
the equality ^A(r,{m e }) = Jr(m e ) holds. 

Let us start by a couple of remarks on the spin foam amplitude. First, the BF partition 
function is topologically invariant, meaning that two triangulations with the same topology will 
yield the same spin foam amplitude. Then, following the analysis of [10, 3], we know that the spin 
foam amplitude is in general ill-defined and that it requires a gauge-fixing in order to give a finite 
amplitude. Roughly, the origin of the divergence is that the S function appearing in the spin foam 
partition function are redundant and, using the topological invariance, we can generally choose a 
smaller set of conditions to impose that the curvature is flat. The completely gauge- fixed amplitude 
for an arbitrary triangulation A corresponds to the smallest triangulation of the same topology and 
compatible with the boundary data (and the graph T being considered as boundary). 
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The main result is that the equality is achieved for any 3d triangulation A topologically equiv- 
alent to the 3-sphere: 

Z A (T,{m e }) =/ r (m e ). (15) 

A general proof in the non-abelian case will be presented in the next section. Let us still explain 
the reason why it works. The reader can also find explicit examples in appendix. 

We start with a graph T embedded in A, and we introduce a framing i.e a surface S C A 
in which T is faithfully embedded. First we consider the case with no non-contractible cycles, 
S ~ S 2 . The triangulation of the surface S directly provides us with a triangulation of 1S3: £3 can 
be realized by taking two copies of the 3-ball glued back together along S2. Then simply choosing 
as triangulation A the initial triangulation for S, the two expressions (10) and (14) simply match. 

The case with cycles is a bit more subtle. Once again, we start with the graph T drawn on 
the surface S, on which we have chosen certain sequences of dual links e to represent the non- 
contractible cycles. Looking at the Feynman evaluation, we associate a variable uc to each cycle. 
To match the spin foam amplitude (14), we associate a face of the 3d triangulated manifold A 
to each cycle: the boundary of this face being the given cycle. More precisely, the cycles of an 
orientable surface S can be split into a set of cycles ai and their dual cycle b{. We construct the 
3d spacetime by taking two copies of S, one for which we add faces corresponding to the cycles en 
and the other for which we turn the cycles h L into faces. The added faces might not be triangles 
at first, but we can triangulate them in the obvious way. Gluing these two filled- up copies of S 
results into a 3d triangulation once again topologically equivalent to the 3-sphere S3. The idea is 
that the interior of S with added faces along one set of cycles, a or b, is topologically a 3-ball. The 
point is that the triangulation resulting from gluing the two copies of S with added faces makes 
the equality between (10) and (14) straightforward. 

Here we have found one particular triangulation Ao for each Feynman diagram T such that 
Zp(Ao) = ir- Nevertheless, due to the topological properties of the spin foam models for BF 
theories, it follows that the same statement is true for any triangulation A topologically equivalent 
to Ao i.e which can be constructed out of Ao by a sequence of Pachner moves (without modifying 
the graph T). 

Let us conclude that this equality between the spin foam amplitude and the Feynman evaluation 
allows us to identify the momentum p e of the particle living on an edge e of the triangulation with 
the holonomy u e = Uf* around the dual face to that edge: this shows us how to encode a particle 
as geometrical data. 

Let us also point out that the Feynman evaluation corresponds to the simplest topology - the 
one of the 3-sphere. The spin foam framework allows to generalize these Feynman graph amplitudes 
to arbitrary topologies. It should be interesting to understand better what effect the non trivial 
topology of the ambient manifold has on Feynman graph evaluation. 

In the next section, we will generalize our framework to the non-abelian context of the Ponzano- 
Regge spinfoam model for 3d quantum gravity, and analyze the Feynman evaluation which results 
from inserting particles in the quantum gravity theory. This will lead us in the following section to 
show that the no-gravity limit of the Ponzano-Regge reproduces the abelian spinfoam model and 
the usual classical Feynman graph evaluation and to identify the Ponzano-Regge amplitudes as 
providing a perturbative expansion in the gravity coupling which is interpreted as QFT amplitudes 
on a non-commutative geometry. 
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3 Particle Insertions in Ponzano-Regge Spinfoam Gravity 



In [3] the general form of Feynman graph amplitude for spinning particles in the Ponzano-Regge 
model has been written. 

In this section we focus on the case of spinless particle (a discussion of the spinning case 
is included in the appendix), we first recall briefly the general construction and then compute 
explicitly the Feynman diagram amplitudes of particles coupled to three dimensional Euclidean 
gravity in a form allowing us to take the "no gravity" limit i.e where the Newton constant Gn — ► 0. 
The main result of this section (see eq. (39, 40)is a explicit computation of the Ponzano-Regge 
amplitude which allow a comparison with the usual Feynman graph amplitude computed in the 
previous section. A similar result has been very recently obtained independently, in the context of 
the Turaev-Viro model, by Barrett et al. [8]. 

We start from a triangulation A of our spacetime M and consider also the dual A*: dual 
vertices, edges and faces correspond respectively to tetrahedra, faces and edges of A. We choose 
our Feynman graph to be embedded in the triangulation A such that edges of V are edges of the 
triangulation. Each edge of T is labelled by an angle 9 € [0, ir] 

9 = Km, k = AttGn (16) 

where Gm is the Newton constant, k is the inverse Planck mass and m is the mass of the particle. 

We choose a group G, here SU(2), and assign group elements g e * to all dual edges e* of the 
triangulation. We constrain the holonomies around dual faces /* ~ e to be flat if there is no particle 
and we constraint it to be in the conjugacy class 9 e if e is an edge of V. More precisely, let us note 
G e = Gf* the product of the group elements around a dual face (or plaquette) /* ~ e: 

Ge = G/.= n 9j: ie '\ 

where ej*(e*) records the orientation of the edge e* in the boundary of the (dual) face /*. The 
amplitude is well defined once we chose a gauge fixing. In order to do so we choose T a maximal 
tree of A \ T and T* a maximal tree of A* [10]. Then the partition function reads: 

Z M {T e )= I [J d 9e* II S(G e )l[A(9 e )6 ee (G e ), (17) 

e*£T* e^Tur eer 

where dg is the normalized Haar measure and 5(g) the corresponding delta function on G and 
A(9) = sin(0). The partition function contains a factor f] e A(0 e ) which can be factor out, this 
factor is important when we consider the no gravity limit. In this section in order to simplify 
notation we will work with the reduced partition function 

ZM{m = mw 

Also, 5g(g) fixes the group element g to be in the conjugacy class 3 labelled by 9. This fixes a 
non-zero deficit angle around the edge e, which corresponds to the geometrical picture of a particle 
in a 3d spacetime of mass AirGm = 9 (e being the trajectory of the particle). More precisely, the 



3 The Cartan subgroup H of SU(2) is the group of diagonal matrices 

h e - 



e w 
e~ u 
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distribution Sg(g) is defined by: 



V/, f dgS e (g)f(g)= [ dx/ixhex' 1 ), (18) 
Jg Jg/h 

where dg and dx are normalized invariant measures. The amplitude (17) does not depend on the 
choice of the triangulation and the gauge fixing trees, but only on the topology of the manifold M 
and the embedding of T in it [4] . 

We can expand the 5 functions in terms of characters 

j 

Se(g) = ^j(he)xM^^x j (0)^x j (he) = S ^-, (19) 

olll (7 

3 

dj = 2j + l being the dimension of the spin j representation. We eventually perform the integration 
over g e * in order to obtain a state sum model 

zm(to) = £ n dj. n km.) n ^ n { f 1 T 2 T 3 1 > (2o) 

{j e }e$r eer e&T t L Jet 4 Jet 5 Jet 6 J 

where the summation is over all edges of A and the product of normalized 6j symbols is over all 
tetrahedra t. For each tetrahedron, the admissible triples of edges, e.g. (je tl ,jet 2 ijet 3 )i correspond 
to faces of this tetrahedra. The factor Sj : o comes from the gauge fixing, it eliminates the sum over 
je, e G T. 

Let us underline the fact that the spin foam amplitudes, like (20), are purely algebraic construc- 
tions built with dimensionless quantities. The gravitational coupling constant Gn does not directly 
enter into the partition function. It appears as a unit used to turn the dimensionless quantities 
(like the angle 0) into physical properties of the matter and particles (like the mass m) . 

Starting from the definition (20) we can explicitly compute these amplitudes. First, in order 
to make contact with usual Feynman graph amplitudes we have to restrict the topology of the 
spacetime to be trivial, so M = S 3 . Also, to warm up, we first consider the case where T is a 
spherical graph which can be embedded in S 2 . Lets start with the example of T being a tetrahedral 
graph embedded in S 2 C 5 3 . We denote by I = 1, ..,4 the vertices of this graph and e = (IJ) the 
edges of this graph. Since the amplitude does not depend on the triangulation we are free to chose 
it. The simplest triangulation of S 3 we can chose in which the graph can be embedded consist of 
two tetrahedra. One of the tetrahedra gives a triangulation of the interior of S 2 the one other gives 
a triangulation of the exterior of S 2 the three sphere is obtained by gluing the two balls together. 
With this graph and triangulation, no gauge fixing is needed, and the corresponding amplitude 
reads 



which correspond to rotations of angle 29 around a given axis (usually the z axis). Every group element is conjugate 
to such an element. The residual action on H — U(l) is given by the Weyl group. It is Z2, since he and h-$ are 
conjugated. 
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The 6j square can be written as a group integral 



U t ;:;F (22> 

where we have introduce the notation J/j = jkl, with /, J, .fT, L all distinct, jjj label the edges of 
the graph dual to r in S 2 . We also denote Ojj = 9kl- Using the previous evaluation and character 
expansion (19) we can perform the summation over spins and obtain 



Zs^e) = f\[dgi J] Sg^gjgJ 1 ). (23) 

J I KJ 

This can be explicitly evaluated [11] 

Z s s(T e ) = Z s s(T e ) n sin^/j) = U-==L_. (24) 
Y<j 2 V det ( cos Oh) 

Lets now consider a general spherical graph T C S 2 . In general T is not a triangulation of S 2 
but we can add edges to it in order to have a regular triangulation of S 2 , we denote it Ta- Once 
this is done we write as before 5 3 as the gluing of two 3-balls, one being the interior of S 2 the other 
the exterior. We can extend the triangulation of S 2 to a triangulation of the interior ball by 
adding interior edges and interior vertices and gluing corresponding tetrahedra 4 . We can choose 
the triangulation of the exterior ball to be the same as the interior ball with reversed orientation. 

Starting from the amplitude (20), we can perform the summation over all edges of the triangu- 
lation that do not belong to S 2 , we obtain 

^(r ) = EII d >n^(^)l^ 3 l r A.i e )l 2 ( 25 ) 

{je} e£T eer 

where the summation is only over boundary edges and (B 3 \T,j e ) is the Ponzano-Regge amplitude 
associated with the 3-ball, with (T, j e ) on its boundary 5 . This amplitude can be understood as the 
physical scalar product between a spin network state and a 'Hartle-Hawking' state associated with 
a manifold with boundary. This amplitude can be easily computed and is given by the evaluation 
of the corresponding spin network. More precisely, let's denote T the graph dual to T in S 2 and Ta 
the graph dual in S 2 to Ta- The edges of Ta are denoted e and are in one to one correspondence 

4 For instance, let's choose an order of the vertices V = 0, 1, • ■ ■ ,n such that vertices / and I + 1 belong to 
neighboring triangles (sharing an edge). We choose as our reference vertex, then we consider the vertex 1. If 01 is 
already an edge of our boundary triangulation we do nothing and consider the vertex 2. If 01 is not an edge of our 
triangulation we add it in the interior of the ball and consider the tetrahedron which consists of this new edge and of 
the two neighbor triangles containing and 1. And then we continue the procedure. This will provide a triangulation 
of the 3-ball which extends the triangulation of S 2 and possesses no interior vertex. 

5 We can extend the definition of Ponzano-Regge amplitude to the case of manifold with boundary dM = E with 
a colored graph drawn on the boundary triangulation to be 

wr,je)=En^n^.oii{ t n t t2 f 3 }• (26) 

Ue} e e ST t ^ JSt 4- Jet 5 JSt e > 

where the summation is only over internal edges. 
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with edges of Ta- Given a coloring j g = j e of the edges of we can consider the spin network 
functional 

*<r Al i.)0te) (27) 

which is a gauge invariant function on with the gauge group acting at vertices. Then 

(B 3 \r A ,j e } = $ {TAh) (l). (28) 
The modulus square of this amplitude can be expressed as an integral 

l^rW 1 )! 2 = / n^n^CftW^J)). (29) 

J v e 

where s(e),t(e) denote the starting and terminal vertices of the oriented edge e. We can now 
perform the summation over jg in (25) in order to get 

Zs^e) = J II d 9v II s eMe)97&) Yl^ma^))- (30) 

" eer e<£r 

We can easily integrate out all the delta functions associated with edges of Ta not in V in order to 
finally obtain 

Zss^e) = [ n ^ n s oMe)97i))> ( 3i ) 

J v£F eer 

which is the desired result. 

We now consider the general case of a graph which is not necessarily spherical but which can be 
embedded in a Riemann surface of genus g, T C As before we can add edges to T in order to 
have a regular triangulation of S g , denoted Ta- Once this is done we write 5 3 as the gluing of two 
handlebodies of genus g: S 3 = H g fe g H*. The meridians of the interior handlebody draw a set of 
a,i(i = 1, • • • , g) cycles on £ ff , the meridians of the exterior handlebody draw a set of bi(i = 1, ■ • ■ ,g) 
cycles on The a cycles intersect transversally the b cycles and together they form a base of 
-?/ 1 (S 9 ). We can extend the triangulation Ta of S g to a triangulation of the interior handlebody 
H g . In order to do so we first have to choose a representative of each meridian as a cycle of edges 
in Ta- Each edge e G Ta belongs to one of the cycle a. L or to none. We introduce an index i(e) 
to keep track of this, where i(e) £ 0, 1, ■ • • , g is equal to if it doesn't belong to any cycle and it 
is equal to i if it belongs to the cycle cij. Each meridian dj is the boundary of a meridian disk Di 
cutting the handlebody. We choose a triangulation A« of Dj which matches the triangulation of a.; 
on its boundary. Once this is done we cut H g along the disks Di and consider H g — UiDi. This is 
a three ball, its boundary contains 2 copies of each meridian disk Di which have to be identified to 
reconstruct H g . Moreover, the triangulation Ta and the triangulation of the meridian disk induces 
a triangulation Ta U Aj U A* of the boundary S 2 of this ball. We can, as previously, extend this 
2-dimensional triangulation to a three dimensional triangulation of the ball. Gluing back the ball 
along Di we obtain a triangulation of H g . For the exterior handlebody H* we do the same except 
that we have to exchange the cycles aj with the cycles b{. In this case we denote i(e) the index 
specifying which b cycle e belongs to. 

Starting from the amplitude (20) we can perform the summation over all edges of the triangu- 
lation that do not belong to we obtain 

Z s ,{T e ) = Y,\{dj e J{x ] Ahe e ){^T,3e\H;){H g \^ v ,j e ) (32) 

{j e } e£T eer 
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where the summation is only over boundary edges and {H g \T,j e ) is the physical scalar product 
between the spin network state |T, j) and the 'Hartle-Hawking' state \H g ). This scalar product can 
be computed, it is given by 

(H 9 \r A ,j e ) = [ flda^ {TA ^(a^) (33) 

where, 3>(r A j^ide) * s * ne s P m network functional, i(e) label which cycle e belongs to, it is if it 
belongs to none, and ao = 1. e(e) = +1 if e as the same orientation than a« and e(e) = — 1 otherwise. 
The proof of this evaluation goes as follows: the idea is to express, in terms of amplitudes, the fact 
that H g is the gluing of a ball along the meridians disks. The edges of the triangulation Ta of H g 
are colored by spins j e , the triangulation A, of Di carries additional edges Sj colored by je r The 
triangulation of the ball Ta U Aj U A* is then colored by spins j e ,jei with the additional condition 
that the spin coloring Aj and A* are the same. The handlebody amplitude is then obtained by 
summing the ball amplitude over all spins 

(H g \T A ,j e ) = ^d^(B 3 \r A UA, l UA*,j e ,j^ (34) 

We know that (5 3 |Ta U Aj U A*,j e ,je t ) is just given by the evaluation of a spin network based on a 
graph dual to T A U Aj U A*. Moreover, this spin network is such that each spin appears twice. 
We can therefore express the evaluation as an integral over group elements associated to the disk's 
vertex v which are in the interior of the -Dj's. 

(5 3 |r A U A, U A*,j e ,je t } = jU dg^ {TAyh) (g He) ) f[ Xjii (st^ 1 ), (35) 

where v{e) is the disk's vertex which belong to e. If e doesn't intersect any meridian disk it is 
understood that gyt e \ = 1. Now the summation over the spins in (34) produces a 5 function for 
every disk edge. One 5 function per disk is eliminated by the gauge fixing. We can then integrate 
out all the group element ga except one per disk which we call aj. This gives us the expected 
formula (33). If we insert this evaluation into (32) we get 

i=i OeWr eer 

This expression is somehow reminiscent of a string theory amplitude 6 . We can now express the 
product as an integral over a product of characters for each edge of T A and then perform the 

6 In the case of string theory the genus g amplitude is written as 

Z= f dmdm^Vi{m)W(m), (37) 

where the integral is over the moduli space of Riemann surface and m is the holomorphic moduli and I labels the 
space of Holomorphic Virasoro conformal blocks. The expression (36) have the same general form if we exchange the 
moduli m, fh by o», bi, the label I by spin label and the Conformal block * by the spin network functional $. It may 
be only an accidental analogy. 
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summation over the j e in order to obtain 

ZMTe) = [ II dg f f flda t db t ]J Se^a^g^) U 6(9^9-^), (38) 
J /cs 9 J i=i eer e ^r 

where we have switched all the notations back to the triangulation so that t e and s e are the two 
faces adjacent to the edge e (they equivalently are the dual vertices v ending the dual edge e). If 
one chooses the cycles a, b to lie entirely in T, we can integrate out the delta functions associated 
with edges not in T and obtain the result we are looking for: 

Z S 3(T e ) = f dg f I f[da t db t n^XS^g))- (39) 
J /cs 9 J i=i eer 

This is the generalization of the Feynman evaluation (10) to the non-abelian case when inserting 
particles in 3d quantum gravity Further assuming that each edge e belongs to a unique cycle of 
T, g , one can simplify the formula : 

Z s s(T e ) = f Udg f f dJJ c Il^e (dteU^g- 1 ) . (40) 

J /c£ 9 J C cycles of S e£T 

4 The no-gravity limit and a perturbative expansion in G 

4.1 The Quantum Field Theory limit of Spin Foams at G — > 

Now we would like to study the "no-gravity" limit of the particle insertion amplitudes of the 
Ponzano-Regge model. That is the limit when we take the Newton constant to zero i.e. Gn — > 0. 
In three spacetime dimensions, the Planck length reads lp ~ HGn while the Planck mass is mp ~ 
1/Gjy- The usual classical limit is taking h — > while keeping Gn 7^ 0. As rap 7^ in this limit, 
we get an effective deformed Poincare algebra [12] and recover the framework of Doubly Special 
Relativity. Here, we consider the alternative limit Gn — ► while h is kept fixed, in this limit, we 
expect to recover from a quantum gravity theory the usual quantum field theory framework: we 
call it the QFT limit. More exactly, we want to recover the usual Feynman diagram evaluations of 
quantum field theory on a flat background as described in the section 2. 

More precisely, the limit is defined as k — > where k is the inverse Planck mass defined in (16) 
as the ratio between the deficit angle and the corresponding mass of a particle. To understand how 
the different quantities gets renormalized in the limit we parametrize g G SU(2) in terms of a Lie 
algebra element In the fundamental representation we have: 

g = e muVi = cos(mK) + i sm(mK)n.a, (41) 

where m = \u\, n is the direction of the rotation and a the Pauli matrices, in the limit, the angle 
6 = rriK goes to so that we are in fact considering perturbations around the identity in SU(2) or 
equivalently the limit in which the group SU(2) goes flat and becomes the abelian group defined 
by its Lie algebra su(2). So it is natural to also call the QFT limit the abelian limit for spin foams. 
We can expand group elements and their products perturbatively in the parameter k: 

( 2l""i| 2 \ -* - 
gi ~ 1 - k —— + iKUi - a + 
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gig 2 ~ 1 + ik (ui + u 2 ) ■ a + — {u\ + u\ + 2u\ ■ u 2 + 2i{u\ A u 2 ) •<?)+... , (42) 

so that the group multiplication is linear at the first order in k. In the k — > limit, the discrete 
representations label j becomes a continuous length parameter I parameterizing the representations 
of 5u(2) ~ M 3 as a group and m is the renormalized mass. The QFT limit is defined as by: 

dj = — , 8 = rriK, k — ► 0, Z,m fixed (43) 
In this limit the integral over the group becomes integral over the Lie algebra 

/ «™ 3 IS < 44 > 

where the normalization factor insures that 

— tf(u) = — [ drum 2 [ d 2 nf(mn), (45) 
2vr vr J J S 2 

with dn the normalized measure on the sphere . Summation over j becomes an integral 

1 r°° 

£ ~ K _o - / dl, (47) 

and we recover the usual classical Hadamard propagator (2) as the abelian limit of the SU(2) 
character. 

, . sindjO 1 sinful , A . 
sm8 k \u\ 

The Ponzano-Regge partition function (17) is given by 

Z M {Te)= [ II d ^ II t(Ge)llS8 e (Ge)A(6 e ). (49) 

Let's first consider the case of a closed manifold without particles and expand the delta functions 
in terms of characters, this gives 

Z m = EIH« / IW ]\xj e {Ge). (50) 

j e e e* e 

Now taking the limit k — ► is straightforward. One uses the fact that the product g e = gf* = 
Hie* ^e* ^ s an aDe h an product at the first order in k then takes the limit of the characters. This 
leads to: 

z - m, z =/n§/„ w. n ^fp - ./», n««- nw. < m > 



7 The normalization comes from the Weyl integration formula 



/ dgf(g) = - f d6sm 2 (6) [ dxfixgx' 1 ) (46) 



with d<j, da; being normalized invariant measures 
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where the vectors u e * are the Lie algebra counterpart of the group elements g e * : 

u e * = m e *n e *, u e = itf* = ti e *- 

We see that we naturally recover the abelian model studied previously. In the limit we have to 
keep track of the factor k. For each edge e not in T, there is a factor 1/k 3 (one coming from ^ •, 
one from dj and one from \j) an d for each dual edge e* not in T* there is a factor k 3 . Overall 
this means that k comes with the power 3%(M), where x(M) = ao — o\ + 02 — 03 is the Euler 
characteristic of A and <7j are the number of i-cells of the triangulation. For a closed three manifold 
x(M) = so no factor k appears in the limit of the partition function. 

A particle insertion of mass 9 e on an edge e, or dual face /*, corresponds to inserting the 
following observable [3] in the partition function (50): 

O e (6 e ) = 8 je0 x Se e (g e ) * A(0) (52) 
where A(6>) = sin 6 and the distribution 5g(g) is defined in (18). Inserting this operator in Z relaxes 
the constraint 5(g) by replacing it by the weaker constraint 5g(g), which gives back the formula 
(17). 

For a particle insertion on the edge e with mass m e , we first observe that the limit of S$(g) is 
8 m (u), which is a distribution satisfying the condition (obtained as the limit of (18)): 

TT"! S m (\u\)f(u) = [ d 2 nf(mn) (53) 

27T Z K° J S 2 

so 5 m (u) = 2 Jm 1 ^\ti\ — m). Note that we have the analog of identity (19) 

f°° sin Im sin I \ u \ ir 

I dl m \u\ =2^ 6i ^- m) - (54) 
The term 5jo in (52) kills the summation J2j djXj(G) which becomes 1/k 3 J dl l sm ^ in the abelian 

limit. So the abelian limit of 5jo is k 3 -jt-- Eventually, the abelian limit of A(#) is Km. This shows 
that the particle insertion in the k — > corresponds to the following observable of the abelian case: 

O e (m e ) = x ^-5(\u e \-m e ). (55) 

Finally, putting everything together, the QFT limit of the Ponzano-Regge amplitude with a particle 
graph {r, m e } is simply equal to the amplitude for particle insertions in the abelian partition 
function and reads: 

Z ( f(T,{m e }) = limK^Z r (9 e ) 



J 



where |er| is the number of edges in V. 

This shows that the QFT limit k — > of the Ponzano-Regge amplitude reproduces the usual 
classical Feynman graph evaluation (with Hadamard propagators), which was previously shown to 
be equal to the amplitudes of the abelian M 3 spin foam model of section 2. The next step is to 
expand the non-abelian spin foam amplitudes beyond the first order in k in order to understand the 
structure of the Feynman graph evaluation in the non-abelian context or, in other words, extract 
the quantum gravity correction to the classical QFT Feynman graph evaluations. 
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4.2 Expansion in G: Non-Commutative Geometry and Star Product 
4.2.1 Feynman evaluation of spherical graphs 

We start with the Ponzano-Regge spin foam amplitude for a particle graph written in section 3. We 
first deal with the case of a spherical graph F Applying the (topology) duality transform described 
in section 2, we can write this amplitude exactly as a Feynman graph evaluation similar to (6). 

More precisely, if T is a spherical graph S 9 ~ S 2 , there is no cycle variables a«, hi to integrate 
over and the spin foam amplitude (39) simplifies down to (31). Then let's denote G e = gt( e )9~^ 
the group element associated to each edge e 6 T, where t(e), s(e) label the two triangles bounding 
e. It is clear that the product of such G e variables around any given vertex v G T is constrained to 
be the identity. And in the end, we obtain the following Feynman graph evaluation: 

l(T,e) = Z s s(T e ) = [ Y\dG e A(G e )8e e (G e ) T[S(G V ), (57) 

with ^ 

G v = l[G?V, 

eDv 

where e„(e) = ± whether v is the target or source vertex of the edge e. Now one can re-express the 
5 function on the group SO(3) as 8 : 

S f G ) = J— f d 3 Xe^ Tr ^\ (58) 
Therefore the non-abelian Feynman evaluation reads: 

= /ilS / T\dGeA{G e )5 0a (G e ) Ue^ T ^\ (59) 

v e v 

If the product G v = Y\ ee g v G e was abelian, we could expand the exponent and reorganize the 
product over vertices into a product over edges just like the usual Feynman graph evaluation (5). 
Nevertheless, we know that the product is actually abelian at the first order in k, so that we can 
do a perturbative expansion in k and compute deviations from the classical Feynman diagram 
expressions. 

Denoting 2inP = Tr(ga) the projection of the group element g on the Lie algebra, the product 
of n group elements g±, ..,g n admits a simple expansion in terms of n: 



in 2 



^PiAPj ■<?+... (60) 



i<j 



8 This formula gives us the delta function on 5*0(3) = SU(2)/Z2, where Z2 is the identification g = —g. In the 
following we restrict to the case where all group elements are SO(3) group elements for simplicity. As pointed out in 
[3] if we want to construct the delta function on SU(2) the correct formula is: 

S(G) = ^-f d 3 Xei Tr< - X3 \l + e(g)), 

where e(g) is the sign of cos 6 with g — cos 6 + i sin 6 — h ■ a. 
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Then we can expand the exponent in the graph amplitude: 



i 

2k 



Tr{X v G v ) = -iX v . ]T [e)P e + K€ v (e) (e') £ ?e * Pe> + ...), (61) 

\e£dv e<e'£dv / 



where X v = iX v .a. So the full Feynman evaluation reads: 

i(r, m ) = /nS / u dG e A ( G ^^ n^-^ 66 ^ 6 ^ 6 ^^^ 6 '^-'^^^- 

v e v 

f e f \ e,e' / 

We also need to expand the measure: 

| ]JdG e A(G e )5 Kme (G e )^G e )^ J n^^lPel 2 - (^^) 2 ) 



in order to express everything in terms of the moment vectors P e and get the full perturbative 
expansion in k. This leads to a simple renormalisation of the mass m — ► sin ^ m , so that the measure 
doesn't contribute at first order in k. If one introduces a source J e for the variables P e , one can 
generate the polynomial terms in P through some derivative with respect to the source J: 



| [— ^ \l + iKX v . V e v (e)e v (e')-^A-^ + ...).I(r,m e ,X v ,J e ) 



j=o 
(63) 

where /(r,?™,^, J) is a normal (abelian) Feynman graph evaluation: 

I(r,m,X v ,J) = jU^ 6 (|P e | 2 - !i^£^ ^ ( p e) JJ e -^e.(x t(e) -x s(e) )-a e .P e _ (64) 

Therefore we see that the Feynman graph evaluation for particles in 3d quantum gravity can be 
expressed as a perturbative expansion in k with operators acting on a classical Feynman graph 
evaluation. 

At first order the operator to insert is a sum of operator acting at the vertices of the Feynman 
diagram which modify how the particles interact. At higher orders the measure will also contribute 
and modify the propagator of the theory. 

One can compute the higher order correction and this perturbative expansion looks quite cum- 
bersome at first sight. However we are now going to see how the full perturbative expansion can be 
simplified if one re-express the previous formulas in terms of a star product, which renders explicit 
the non-commutative structure of the theory and scattering amplitudes. 

4.2.2 Star Product and Non-Commutative Field Theory 

Starting from the non-abelian Feynman evaluation (59), one would like to split the exponentials 
ex.p(^Tr(X v fl e G e )) in order to re-organize all the exponentials and write the evaluation as a 
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standard Feynman amplitude (as in section 2). It is thus natural to introduce a notion of ^-product 
on functions of X G su(2) ~ B 3 denned through: 

e ±Tr(X 9l ) ^ e £Tr(Xg 2 ) = e £Tr(X 5l92 )_ (g5) 

This leads to introduce the following Fourier transform between the Lie algebra su(2) and the Lie 
group SO(3): 

f(X) = J dge^ Tr ^f(g). (66) 

This ^-product is non-commutative but still associative. Moreover it can be seen as equivalent to 
the convolution product on the group: 

(^^)=0o G ^, with 0o G ^( 5 ) = J ' dh^gh- 1 )^). (67) 

We can also define the star convolution product on R 3 which is dual to the usual product of 
functions on the group. 

(J^J)=00, with <j)o^(X)= f d 3 X<fi(X — Y) * ip(Y). (68) 

JR3 y 

A natural set of functions to study are the Fourier transforms fj of the character Xjid)- It is eas Y 
to derive that: 

fj*fk = ^fj, (69) 
j+k 

i=\j-k\ 

More generally, if we denote D ] ab (g) the matrix elements in the spin j representation, we can define: 
we can check that: 

fL * fed = ~^~^bcfid 

= E C lceCti}f ef , (71) 
l,e,f 

where we have defined f^J{g) = D^ d (—g~ l ) and the Clebsh-Gordan coefficient C jkl of the recoupling 
theory for representations of SU(2). This shows that the f/s project onto the j representation for 
the ^-product: 



f . fk _ °Jk f j 
J]~* Jed — ^ Jed' 



In particular, f is a projector: 

/o*/^ = <5 ofe /o, (72) 
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and plays the role of the usual S(x) distribution. In order to understand the behavior of the 
functions fj, it is useful to look at the link between this algebra/group Fourier transform and 
the standard Fourier transform on M 3 . We can express SO(3) group elements in terms of their 
projections on the Lie algebra 



g = ^1- k 2 \P\ 2 + inP.a (73) 

A function / on SO (3) can be equivalently seen as a function on the 3-ball B 3 = {\P\ < k' 1 } 9 . 
We call this space of function B K (M. 3 ). This is our Fourier space: by construction no momenta in 
this space takes a value larger than re -1 . Our group Fourier transform is related to the usual fourier 
transform by: 

/(*) = £/ / ^ e-*- p f{P). (74) 



The functions fj can be seen to be related to the Bessel functions of the first kind. The important 
feature is that /o has a non-zero width. Therefore, as /q plays the role of 5(X) for the ^-product, 
this means we have access only to a finite resolution: the ^-product tells us of a minimal length 
scale (maximal resolution) on the spacetime (X sector). This is simply due to the fact that the 
momentum space (P sector) is bounded. 

In order to describe the functional space dual to ,6 K (IR 3 ), we introduce the following kernel: 



This is a projector as we can check that J d 3 X G(X,Y)G(Y, Z) = G(X,Z). It is now straightfor- 
ward to show that if f(X) is the group Fourier transform of a function / then G(f) = f. Moreover 
any function in the image of G has a Fourier transform with support in B 3 , therefore the image of 
L 2 (M 3 ) by this projector, denoted B K (M. 3 ), is isomorphic to B K (M. 3 ) by the group Fourier transform 
(74). This Fourier transform is in fact an isometry if we provide B K (M. 3 ) and B K (M. 3 ) with the 
following norms: 

2 r d 3 x 



= 4 /. 1 ^ l^)! 2 - (76) 
With the tool of this star-product, the Feynman evaluation (59) reads: 

j ( r ^) = /nS / n^ A (^)WGy n ( e * e ^ Tr(XuG " (e)) ) • (77) 



9 If f is continuous at the identity, it satisfies the additional condition f(n 1 n) is independent of n a unit vector, 
we will not impose this condition in general 
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When written in terms of momenta 10 

m n - (n / n S /n^%p * - n (*/*<•*• 

2 (81) 

So far the propagator which enters this amplitude is the Hadamard propagator 5 (\P e \ 2 — sm2 JT eK ^ • 
We can write this propagator as a proper time integral: 

V k 2 y Loo 2vr 

If one restricts the integral to be over positive proper time T € K + , one obtains the usual Feynman 
propagator with a renormalized mass. Using this Feynman propagator, the spin foam amplitude 
reads 

*™ -nfi?)/ng/n* , Ap .Jw* ? (* — *) ■ <«> 



with P e = P(g e ). 

In the next section, we will show that these amplitudes are truly the Feynman graph evaluations 
of a non-commutative field theory for a /t-deformed Poincare group. Such field theory then acquires 
the interpretation of an effective theory for matter propagating in the 3d quantum geometry. 

4.2.3 Non-spherical Graphs and non-trivial Braiding 

Beyond the case of spherical graphs, dealing with graphs with an surface embedding of non-trivial 
topology g / is more subtle. It will lead us to the notion of braided Feynman diagrams for a 
non-commutative field theory. 

Starting with the graph amplitude (39) and defining the edge group elements H e = g^al^ g~^b £ e ^ 
it is obvious that the product of such H e variables around a vertex v 6 T will not generally be 
the identity. First, let us notice an ambiguity in the definition of the group element associated 
to an edge. Indeed the distribution Sg is invariant under conjugation Sg(-) = 5g(k ■ k^ 1 ), so that 
we can in general choose some arbitrary variables k e and use the group elements G e = k e H e k~ l . 
To fix this ambiguity and build rigorously the group variables G e , one uses a methods similar to 
the gauge-fixing techniques used in [10, 3, 13]. To follow the procedure, it is easier to work on 
the dual graph to T. One chooses a dual vertex vq of reference (an initial face /o) and a maximal 
tree T on the dual graph. For any dual vertex v, there exists a unique path V(vq — > v) running 
from vq to v along the tree T. Then one can define the ordered product ky of group elements g^ 

10 The normalizations relating g and P are given by 

(78) 



^ = 2^^ S {^ S3-)' ^ 

A(kw) = sin Km. (80) 
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Figure 1: Case of a crossing: we lift the crossing by embedding it into a handle, we triangulate 
the corresponding punctured torus and we derive the braiding rule in the non-abelian Feynman 
evaluation. 

associated to the dual vertices along the path V{vq — > v). One should moreover include all the 
cycles ctj, bi that the path crosses: if the dual edge between the consecutive dual vertices w and w' 
intersects a cycle (i.e if the corresponding edge belongs to the cycle), then one includes the group 
element corresponding to that cycle in the ordered product k^. Then one would define the edge 
group element G e = k t ^H e k~}y which is simply the ordered product of face group elements gf 
and cycle variables around the edge e but starting at a fixed reference face fo- Nevertheless, the 
(ordered) product of the G e 's around all vertices is still not trivial and the situation is more subtle 
than a simple choice of origin on the dual complex. 

The moot point is at the level of crossings: when two edges of the graph cross each other (on 
a 2d projection of the graph or more precisely when embedding the graph on the sphere), we have 
a non-trivial braiding of the edge group elements, which forces us to associate two group elements 
to each edge, one at the source vertex and one at the target vertex. 

First T being a non-spherical graph means that we cannot draw it on a sphere without crossings. 
In order to compute the amplitude we need to choose a surface on which this graph is drawn without 
crossing. There is an arbitrariness on the choice of this surface of course. Our amplitude being 
topological, it doesn't matter which surface we choose and we are free to choose the most convenient 
one. The surface which we use consists of adding one handle to the sphere for each crossing: the 
upper edge goes through the handle whose meridian is a cycle a and the lower edge goes below the 
arch of the crossing as shown in fig.l. This means that starting from a projection of our diagram 
on the sphere we cut a small disk around each crossing of the projection and glue back a punctured 
torus (with one hole). Topological invariance insures that the amplitude doesn't depend on the 
choice of projection of our graph. 

Then let us focus on one crossing. We choose the simplest cell decomposition of the punctured 
torus as on fig.l. There is 4 faces to which we assign group elements g±, ..,34. In this diagram there 
are 8 edges. Each edge lead to a potential constraint. 4 of them are mass constraints, only two 
(one per edge) are relevant. Among the other 4 constraints, one can be gauged away (Lets say the 
one associated with the edge bounding the face 2 and 3). Then the spin amplitude (39) for the 
concerned faces and edges around the crossing reads: 

<5<?i {91 9a 1 ) $e 2 {gi9 2 1 ) 8 (92 9\ 1 b) 5 (g^ 1 b)5(g 4 ag^ 1 ) 

Integrating over the cycle group elements a, b leaves a single constraint ^{gig^ 1 g^g 4 l ) plus the mass 
constraints. We then introduce the group elements associated to the two edges at each end: 

G S A = g^A = 9392 \G S B = gig 2 \G l B = g A g^ 1 . 
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Figure 2: Feynman rules for particles propagation in the Ponzano-Regge model. 

They satisfy the constraints G s A G s B (G t A )^ 1 (G B )~ l = 1. In terms of these edge variables, the spin 
foam amplitude reads: 

<5 01 (Gi)^ 2 (G|)«5(GiG^(G t A )- 1 (G t s )- 1 ) ( 5(G^(G t s )- 1 ). (83) 

The important thing to notice is that we do not have a constraint d(G A (G t A )~ 1 ). 

We can now fully specify our Feynman rules. We first assign two group elements to each edge 
G s r e \, G t r e \. For each vertex we put a weight ^dleDu ^s(e))> f° r each edge with no crossing we put a 
weight 5(? e (G s ( e ))(5(G s ( e )(Gi( e )) _1 ) and for each crossing we add <5(G s ( e /)G s(e )(Gi( e /)) _1 (Gi( e )) _1 )(5(G s ( e )(G t ( e )) _1 ) 
where e' over-crosses e. This is summarized in fig. 2. We show in the next section that these Feyn- 
man rules can be derived exactly from a non-commutative field theory based on the K-deformed 
Poincare group. 

These Feynman rules can be understood as a Reshetikhin-Turaev evaluation of a colored graph 
constructed with the quantum group T>(SU (2)) (which is a K-deformation of the Poincare group) [4]. 
They involve a non-trivial braiding factor for each crossing. We can do the same trick as previously 
and promote the Hadamard amplitudes to Feynman amplitudes by taking only the integration over 
positive proper time in the propagator and we call the resulting amplitudes Xp(r, 6). 

5 A Non-Commutative Braided Field theory 

5.1 Non-Commutative Field Theory as Effective Quantum Gravity 

We now introduce a non-commutative field theory based on the previous ^-product, and we show it 
leads to a sum over graphs of the particle insertion amplitudes obtained from the spin foam model. 
It has the natural interpretation of an effective field theory describing the dynamics of matter in 
quantum gravity after integration of the gravitational degrees of freedom. 
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Let us now consider the graph evaluations 2f(T, 6) given by the equation (82). We consider the 
case where we have particle of only one species so all masses are equal, m e = m. Having different 
masses would only require to introduce more fields and would not modify the overall picture in any 
way. We define a normalized amplitude 11 2f(T,8) = ( co ^1 m ) ^2f(T, 0) where |er| is the number 
of edges of V. Let us now consider the sum over planar and trivalent graphs: 

-^l F (T^m) (84) 

Tplanar , trivalent 

over trivalent planar diagrams (we will include non-planar diagrams in the next section) where A is 
a coupling constant, \vp\ is the number of vertices of V and Sr is the symmetry factor of the graph. 
Remarkably, this sum can be obtained from the perturbative expansion of a non-commutative field 
theory given explicitly by: 



S =8^ ld6x 



\(di<j> * di<j>) (xf - \ {4> * 0) (x) 2 + A (0 * * 0) ( X ) 



(85) 



where the field 4> is in B K (M. 3 ). Its moment has support in the ball of radius re l . We can write this 
action in momentum space 

S(4>) = \J " dg (p\g) - S -^P^j Mtig- 1 ) + dg 1 dg 2 dg 3 5(g 1 g 2 g 3 )^gS(g 2 )^(g 3 ), (86) 

from which it is straightforward to read the Feynman rules and show our statement. 

It is interesting to write down the interaction term in terms of the momenta P(g). The momen- 
tum addition rule becomes non-linear, in order to preserve the condition that momenta is bounded, 
explicitly: 



Therefore the momentum conservation at the interaction vertex reads: 

Pi © P 2 © ^3 = = Pi + P 2 + P 3 - k(Pl A P 2 + P 2 A P 3 + Pi A P 3 ) + 0(re 2 ). (88) 

From this identity, it appears that the momenta is non-linearly conserved, and the non-conservation 
is stronger when the momenta are non-collinear. The natural interpretation is that part of the 
energy involved in the collision process is absorbed by the gravitational field, this effect prevents 
any energy involved in a collision process to be larger than the Planck energy. This phenomena is 
simply telling us that when we have a high momentum transfer involve in a particle process, one 
can no longer ignore gravitational effects which are going to modify how the energy is transferred. 

The non-commutative field theory action is symmetric under a re-deformed action of the Poincare 
group. If we denote by A the generators of Lorentz transformations and by the generators of 
translations, it appears that the action of these generators on one-particle states is undeformed: 



A = (^(KgAr 1 ) = (p(A ■ P(g)), 
Ta-fcg) = e iP ^- 3 4>(g). (90) 



11 From the spin foam point of view each amplitude correspond to a manifold of different topology. The spin foam 
approach does not specify how one should relatively weight amplitude of different topology and we have to make a 
natural choice. 
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Since there seems to be a certain amount of confusion in the literature on this subject, we would 
like to emphasize that it is impossible to deform non-trivially the action of the Poincare group 
on one-particle states. Indeed, it is well known that the cohomology group of the Poincare group 
is trivial: any deformation of the Poincare group which is connected to the identity (i.e. which 
depend continuously on a parameter k) can be undone by a non-linear redefinition of the algebra 
generators. More precisely, let us start from any /i-deformation (which preserve associativity or 
Jacobi identity) of the Poincare algebra, that is we have Poincare generators J^ v , P^ with deformed 
commutation relations such that when k = we recover the usual algebra. Then we can always 
choose new generators J^ v , P^ which are non- linear functions the original generators and k such 
that the new generators satisfy the undeformed Poincare algebra. This choice of generators is 
referred as the classical basis. This is nicely exemplified in [5] for the 3+1 n- Poincare algebra and 
in [12] for the 2+1 re-Poincare algebra. The open problem in this context is whether physics does 
depend on the choice of basis or not and if so which one is preferred. It seems that the choice of 
basis matters, at first sight, since each basis correspond to a choice of what should be used as an 
operational notion of energy [14], and for instance the dispersion relation is dependent on the basis 
choice. Whether this is a true physical dependance or not is still however a matter of debate [15]. 

It has already been noticed that particles in 3d quantum gravity should be described by a 
Deformed/Doubly Special Relativity (DSR) [3, 12]. What we learn from the present analysis is 
that 2+1 gravity chooses for us a preferred basis which is a classical basis (this fact was already 
noticed in [12]). Note however that if J^^P^ form a classical basis, then any redefinition of the 
form J^y = J^u^P^ = P m /(k|P|), with an arbitrary function /, still yields a classical basis. So 
2+1 gravity chooses one particular classical basis 12 . This chosen one, singled out by gravity, is 
consistent with the DSR principle since it is not possible in this basis to have an energy which 
exceed the Planck energy. 

In a classical basis, the non-trivial deformation of the Poincare group appears at the level of 
multi-particle states. Indeed, the action of the Poincare group on two-particle state is given by 



The action of the Lorentz generator is undeformed however the action of the translations is 
modified in a non-trivial fashion reflecting the fact that the spacetime has become non-commutative. 
With these rules it is clear that our field action is k- Poincare invariant. In algebraic language this 
means that the symmetry algebra is promoted to a non-cocommutative Hopf algebra. In the 
classical basis, the only non-trivial coproduct is 14 



This shows that the effective theory describing the dynamics of second quantized particles in the 
presence of gravity can be effectively described, after integration over the gravity degrees of freedom 
as an explicit non-commutative field theory which respects DSR principles. 

12 If we denote p the standard classical basis and P the one appearing in our context the relation is given by 

pa _ sin«[p| a 
K \p\ V 

This is different from the K-deformation of the 3+1 Poincare algebra. 
14 This is exactly the Hopf algebra structure of the ft-deformation of the 2+1 Poincare algebra, or equivalently of the 
quantum double 2?(SU(2)), which appears in the quantization of 2+1 gravity [4, 17]. The infinitesimal deformation 



A-<KPi)<KP 2 ) = <KA-Pi)<MA-P 2 ), 



(91) 
(92) 



A(Pi) = ^l-P 2 ®Pi + Pi® Vl-P 2 ~ e ijk P j ® P k - 



(93) 
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In fact we will now see by looking at the non-planar graphs that the proper effective description 
is in term of a braided non-commutative field theory which takes into account the non-trivial 
statistics induced by gravity. Also, the Euclidean gravity amplitudes which we have discussed so 
far are expressed in terms of Hadamard propagator and we had to extract the Feynman propagator 
from them, in order to make contact with non-commutative field theory. We will argue in the 
next section that the Lorentzian gravity spin foam amplitudes should have a formulation in which 
gravity amplitudes are given in terms of causal propagators. 

5.2 Non-Commutativity and Braided Feynman Diagrams 

The question is now whether the non-commutative field theory introduced above reproduces the 
quantum gravity amplitude in the non-planar case. The computation of Feynman amplitudes in 
the planar and non-planar case differs in the sense that in the non-planar case we have to commute 
the Fourier modes of the field before doing any Wick contraction. For instance if we compute using 
Wick theorem the following vacuum expectation value, we get 

<to)£(P 2 )fe)feO> = (HPMP2)) G>{Pz)4>{P±)) + (HPi)HPa)) <fe)fe)> 

+ {<f>(P 1 )cf>(P 3 )){ ( f>(P 2 )<f>(P 4 )). (94) 

If we draw the corresponding Feynman diagrams, putting the momenta ordered on a line we see 
that the first two contractions are planar whereas the third one contain a crossing, the crossing is 
due to the fact that we have to exchange 4>{P 2 ) and (j>{P 3 ) before making the Wick contraction. 
In order to compute the non-planar Feynman diagrams we therefore have to specify the rules for 
commuting the modes <ft(P), that is specify the statistics of our particles. In standard commutative 
local field theory, we know that only two statistics are usually possible bosons or fermions (except in 
2 + 1 dimensions). What about non-commutative field theory? It seems that this is terra incognita 
since the usual spin statistics theorem cannot be applied (except for some particular examples in 
usual non-commutative field theory with no space time non commutativity [16]). We now argue 
that once we have fixed the star product and the duality between space and time there is a natural 
way to specify the statistics of our field. Let us look at the product of two identical fields: 

<t>*<t>{x) = y^i^ 2 e^ tr(a:si92) ^i)0(52), (95) 

We can move 4>(g 2 ) to the left by making the following change of variables g\ — > g 2 and g 2 — > g 2 1 gig2i 
the star product reads 

</>*<f>(x) = j dg 1 dg 2 e^ t ^^(g 2 )0(g 2 1 g 1 g 2 ), (96) 

This suggest that the proper way to read the statistics of our non commutative field is to assume 
that they satisfy the commutation relation: 

4>(9i)4>(92) = 4>{92)^{g 2 1 gig2) (97) 

of the algebra is simply given by: 5(Ji) = 0, 5(Pi) = tijkP'' ® P k ■ The K-deformation of the Poincare algebra in 3+1 
dimensions is actually different and usually involves singling out the time direction and also deforming the Lorentz 
generators. 
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In our case, we can check that this commutation relation is exactly the one coming from the braiding 
of two particles coupled to quantum gravity. This braiding was computed in the spin foam model 
[3] and is encoded into a braiding matrix 

R ■ H91M92) = ktoWtfgm). (98) 

This is the R matrix of the K-deformation of the Poincare group [18]. We see that the non-trivial 
statistics imposed by the study of our non-commutative field theory is related to the braiding of 
particles in 3 spacetime dimensions 15 Such field theory with non-trivial braided statistics are usually 
simply called braided non-commutative field theory and they were first introduced in [19]. 

If one uses the non-trivial statistic (97), it becomes an easy exercise to check that the parti- 
tion function of our non-commutative field theory reproduces the sum over all quantum gravity 
amplitudes with insertions of Feynman propagators. 



6 On the classical limit of the Turaev-Viro model 

6.1 The many limits of the Turaev-Viro model 

We have dealt so far with the case of a non zero cosmological constant A. When A / the quantum 
gravity amplitudes are known to be given by the Turaev-Viro model instead of the Ponzano-Regge 
model. We want to show in this section that the zero gravity limit of Turaev-Viro reproduces com- 
putation of Feynman diagram with insertion of the Hadamard propagator on S 3 . We also present 
a modification of the Turaev-Viro model which reproduces, in the no gravity limit, Hadamard- 
Feynman diagrams on the hyperbolic space 

If we consider the case of a non-zero cosmological constant we have three dimensionfull constants 
at our disposal,the planck constant h, the Newton constant k and the cosmological scale L = 1/ y/X. 
We have two lengths at our disposal: the maximal cosmological length and a minimal planck scale 
lp = Hk (keep in mind that we have set the speed of light c to 1). We can therefore introduce a 
dimensionless ratio 

" = TT- (99) 

k is then quantized and labels the quantum group deformation parameter q = exp ^ (see e.g [20]). 
The two length scales also give rise to two mass scales. The physical implication of these mass and 
length scales imply that any change of length or mass Al, Am is bounded from above and below 

Hk< Al < ttL, (100) 

% < Am <-. (101) 
L K 

The bounds on Al imply that I is discrete in Planck unit and bounded from above and similarly the 
mass is discrete in cosmological units and bounded from above. This is indeed what is realized in 

15 It is often believe that this is true only because we are studying three-dimensional field theory and this effect 
should disappear in higher dimension because it is impossible to have a non-trivial statistics in 3 + f dimensions. This 
relies on the fact that in a classical spacetime the homotopy group of the configuration space of n particles on the 
sphere is the pure symmetric group. This belief amounts to suppose that such a theorem still holds true in the context 
of a non-commutative space time. This is far from being obvious and there is no evidence supporting this hypothesis. 
Moreover anybody who has dealt with non-commutative space-times knows that in a deep non-commutative regime, 
the notion of dimension is not sharply denned, and our argument tends to be insensitive to dimensionality. 
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the Turaev-Viro model. It defines a state sum model where the spins j representing the geometrical 
information label representations of a quantum group SU (2) q and the summation over j is bounded 
from above dj = 2j + 1 < k. Then the physical length is related to the spin by 

I = hndp (102) 

The insertion of spinless particles in this model have been described by Barrett in [7] and are 
analogous to the insertion of particle in the Ponzano-Regge model. Instead of inserting A(6)xj(hg) 
in the state sum measure along the edge of the Feynman graph, one puts Sj a where a is a quantum 
group representation label (so d a < k). The spin a is related to the physical mass by 

Ma MOQ\ 

m = — — (103) 
Li 

and S is given by 

sin j;djd a 

S ja = *^-A 104 

3a sinf K J 

Sjo is the quantum dimension. In this model both lengths and masses are discretized and satisfy 
the bounds (100). The amplitude is 

z qM $ ) = n^y: n s >o n n { - et 1 t ta f 3 }' (io5) 

{je}e<£T e£T t ^ Jet t Jet 5 Jet 6 >q 

where we use the quantum group 6j symbol. N is a normalization parameter given by 

"-Xt-t^j- (106> 

\v\ is the number of vertices of the triangulation. One can show that this definition of the Turaev- 
Viro amplitude is equivalent to the gauge fixed definition. 

Starting from the Turaev-Viro model, we can take several limits of the fundamental constants 
while keeping physical length and mass fixed: 

• The zero cosmological constant limit L — > oo: 

In this limit L — > oo the upper bound on length and the lower bound on mass disappear, 
we then recover Ponzano-Regge model studied previously where lengths are discrete and 
unbounded, and masses continuous and bounded. 

• The no gravity limit k — > 0: 

This is the limit we are interested in. In this limit, lengths become continuous and bounded 
while masses turn out discrete and unbounded. 

• The semi-classical limit h — ► 0: 

In this limit both the mass and the lengths are continuous and bounded. It would be inter- 
esting to understand better the resulting state sum model. 

Note that in all these limits the deformation parameter becomes trivial q — ► 1. But the simple 
statement q — > 1 is not enough to specify which limit we are considering. We also need to define 
how the basic quantities like masses and lengths are rescaled in the limit. More precisely, the limits 
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are defined as the same mathematical limit q — > 1, however, physically, we need to specify the 
units for the physical quantities such as the mass and the length: different choices of units lead to 
different effective theories. 

In the following we work in units h = 1 and we consider the no-gravity limit of Turaev-Viro. 
Let x,y be two points on a 3-sphere of radius L and denote by l(x,y) the distance between them. 
The Hadamard propagator on a sphere of radius L is given by 

„ , . smml(x,v) ( * 9 1 \ „ , 
G ^ ) = ^f- (A + ^-^jG^O, (107) 

Where A is the Laplacian on the sphere. If we denote by x, y the corresponding points on the unit 
sphere SU(2) and by a the half integer such that d a = rriL we see that the Hadamard propagator 
is just the character G m (x,y) = Xa{xy~ V )- The evaluation of Feynman diagram with insertion of 
the Hadamard propagator of S 3 is given by 

Z(T,a) = f HdxvUxaMie)^) (108) 

The integrals over x v can be easily performed, and we are left with the evaluation of the spin 
network functional at the identity 

Z(r,a) = |$ (rA , ae) (l)| 2 . (109) 

We claim that this is the no gravity limit of the Turaev-Viro amplitude. For instance if we consider 
the tetrahedral graph as an example we have 

*cr..)=/n*,n*.w)={£ z 2f- (110> 

One can check that this is actually the limit k — ► of the Turaev-Viro amplitude: 

Z q (r,a) = -LTT\S jljai A JU m i23l 2 = |ai2 a 13 a 14 | 2 _ 
q N 3 f^f-\ 3ljaiJ \ j 12 ji 3 314 J I a 34 a 24 a 23 J v ; 

In general, it has been proven by Barrett [8] that Z q (F,a) is given by the square of the quantum 
group evaluation of the colored graph (r, a). The no gravity limit of this evaluation therefore 
reproduces the Feynman graph evaluation (108). 

We are now going to show that the Turaev-Viro amplitude (105) can be written as the Feynman 
graph evaluation of a non-commutative braided quantum field theory. The classical propagator is 
a function on 5 3 x 5 3 . In a non commutative field theory 5 3 becomes a non commutative (fuzzy or 
q-deformed) space and functions on 5 3 becomes operators. If we pick an axis in 5 3 and denote R 
the distance from the north pole of 5 3 along this axis, we can describe 5 3 as a stack of two spheres 
of radius Lsin(R/L). Our fuzzy S 3 is then a stack of fuzzy two spheres. It is well known [21] that 
the space of function on a fuzzy two sphere of a given radius is given by f q {S 2 ) = End(V^), where 
Vj is here the quantum group representation space of spin /. The 'radius' of this q-deformed sphere 
is given by the quantum dimension dim q (Vj) = S1 °^y^ . The space of function on the q-deformed 
S 3 is then given as a stack of fuzzy spheres 

(fc-2)/2 

F q (S 3 )^ End(Vj). (112) 

7=0 
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i> J 



= G a (i,j) 




= R(i,j,k,l) 



Figure 3: Feynman rules in the Turaev-Viro model: a Feynman graph evaluation is given by the 
corresponding Reshetikhin evaluation of the diagram. 

A function is thus given by a collection of operators /(/) 6 End(V/) and the normalized integral is 
replaced by a q-trace 



The propagator G a G F q (S 3 )®F q (S 3 ) becomes a collection of operators G a (I, J) G End(V/®Vj) 
which are given by the Reshetikhin evaluation of the diagram. The Feynman rules are listed in 
figure 3. For each propagator we insert G a (I, J), at each vertex we perform the trace and sum over 
I and for each crossing we insert a R- matrix. This last factor means that, unlike the usual view on 
non-commutative field theory, the braiding of the graph matters for the Feynman graph evaluation 
and we are dealing with a braided non-commutative field theory [19]. This expresses the fact that 
non-commutativity is naturally associated with a change of statistics of the quantum particles. 

It is easy to see by performing the sum over / that our Feynman rules gives the evaluation 



where (r, a) q is the Reshetikhin- Turaev evaluation of the graph T colored by a e . It is less trivial 
to show that this is equal to the definition (105), this is proven in [8]. 

6.2 A hyperbolic quantum state sum 

We have so far considered Feynman graph evaluation of particle propagating in an Euclidean space 
having a positive or zero cosmological constant. For completeness we now consider the case where 
the cosmological constant is negative and the particles propagates in the hyperbolic space H 3 . 

We claim that the corresponding state sum model is a modification of Turaev-Viro where the 
deformation parameter is taken to be real q = exp( — ?), no quantization condition holds on k. 
The representations of U q (SU(2)), q G M, are labelled by unbounded spins j = 0, ..,oo. The state 
sum model is (up to the normalization factor, the range of summation over j and the gauge fixing) 




(113) 



Z q (T,a) = \(T,a) q 



(114) 
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similar to (105) 



;wi»=£ n dim^n^n^oiil - eti f 2 f s 1 - (n5) 

{je} e^Tur eer eST t L *5 Jet 6 j q 

where T is a maximal tree of A\T, dim,j(j) = s ' s ^^/fc)^ ^ s ^ e quantum dimension, p is a continuous 
angle € [0,k] and 

s*-^- die) 



We now have the remarkable property that the no-gravity limit of this amplitude is given by the 
Feynman graph evaluation of Hadamard propagators in H 3 . The Hadamard propagator on H 3 is 
given by 

sin pl(x,y) 

Up[X,V) ~ psmhl(x, y y [ ' 

where l(x,y) is the hyperbolic distance. The Feynman graph evaluation on H 3 is given by 

Z(T,a)= / Y[dx v Y[G Pe (x t{e) ,x s{e) ). (118) 

v e 

This amplitude is an evaluation of a colored graph of SL(2, C). p label simple unitary representation 
of SL(2,C). The integration at each vertex produces a simple intertwiner of SL(2,C) according 
to the general picture described in [22], and the recombination of these intertwiners produces 
the usual group evaluation of colored graph. This evaluation is exactly the so called relativistic 
evaluation used as a building block for Lorentzian gravity amplitude [23]. For a tetrahedral graph 
this amplitude gives the 6j symbol of SL(2, C) 

A detailed proof of this claim follows from the results of [24, 25]. One can show that the state 
sum (115) is an evaluation based on the non-compact quantum group U q (SL(2, C). Namely (p, 0) 
labels the unitary simple representations of ?7g(SL(2,C). Given such a labelling and a colored 
graph (r,p e ), we consider the U g (SL(2, C) Reshetikhin-Turaev evaluation |(T, p e )\, the intertwiner 
between the simple representations at a vertex is constructed in [24]. The no-gravity limit of this 
evaluation is then given by the evaluation of a graph colored by simple representation of the group 
SL(2,C). This proves the claim once we have check that the normalization agree (the interested 
reader can find more details in [25]). 



7 Causality: Hadamard function vs Feynman propagator 

7.1 A brief review of Particle Propagators and Feynman Amplitudes 

If we would like to recover the (scattering) amplitudes defined by the Feynman diagrams of the 
field theory, we would like to attach Feynman propagators to links of the graph T instead of the 
Hadamard Green function. 

Let us first review the definition of the different Green functions. The basic building blocks are 
the positive and negative Wightman functions. There are solutions to the Klein-Gordon equation 
defined as: 

G ± (x\y) = J d 3 p9( Po )5(p 2 - m 2 )e^P< x -y\ (119) 
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and expressed in terms of correlations of the scalar field 4> as: 



G + (x\y) = <0|<K*)<Ky)|0>, G-(x\y) = {0\<j>(y)<f>(x)\0) = G+{y\x) = G+{x\y). (120) 
The Hadamard function is also a solution of the Klein-Gordon equation. Expressed as: 

G(x\y) = J^yj d 3 P 5(p 2 - m 2 )e-^ x -y\ (121) 

it is the sum G = G + + G~ . The Feynman propagator Gf satisfies the equation (D x + m 2 )GF = 
— 5^\x — y) and can expressed as: 

%G F {x\y) = 9(x° - y°)G + (x\y) + 8(y° - x°)G-(x\y), (122) 

iG F {x\y) = <O|T0(x)0(y)|O). (123) 
The proper time expression for G and Gf are: 

G(x\y) = —^ dT / d 3p e -i>(*-»)-r(p a -m 2 )] j (124) 
(27r)- 3 J 

-I r+oo r 

G F (x\y) = —iyz — / / rf3 pe -i|p-(*-»)-T(p -m +ia)] > (125) 

1 /* e-ip-ix-y) 



d A p 



(27r) 3 7 p 2 — m 2 + ia 

where a > is just a regularization. 

Finally, one can also introduce the causal Green function Gc , which is a solution of the Klein- 
Gordon equation. Its integral expression is: 

2 _ ™ 2 V=Fip-(z-j/) 



Gc(z|y) = j d 3 pe(p )5(p 2 - m 2 )e 



and correspond to the correlation: 

iG c {x\y) = (O|[0(x)0(y)]|O) = G + (x\y) - G~{x\y). 

Now let us consider a Feynman diagram i.e a graph T with Feynman propagators living on its 
edges. The evaluation of the diagram is the amplitude: 

I T = f [] dx v J] <?£ (x t{e) - x s(e) ). (126) 

For a given set of positions {x v } and thus a particular (causal) ordering of the vertices, the product 
of Feynman propagators is simply a product of Wightman functions: 

II G mM(e) " f s(e)) = II G ™M(e) ~ %s(e)), (127) 



32 



where the e e 's record the ordering of the vertices i.e the future/past orientation of each edge. This 
way, it may seem natural to introduce Wightman diagrams, with a causal structure {e e } imposed 
independently from the true ordering of the x's, whose evaluation will be: 



= / ]J dSv II G 'k ( f * (e) - ^(e) ) 
•* v e 

]Jd 3 p e [ nd 3 x^(e e p°)<5(p2- m 2) e -^-fe( e )-^(e)) 



/ 



Hd 3 p e 9(e e p e )S{p 2 e -m 2 e )l[s[ £ p e - Y, Pe). (128) 

e\v=t(e) e\v=s(e) 



Because of the constraint imposing momentum conservation at the vertices, one can solve these 
constraints as explained previously and pass to the dual formulation of these amplitudes: they can 
be easily reproduced by a spin foam model. 

In order to do the same with the Feynman amplitude, one should get rid of the explicit use of 
the causal ordering. The solution is to use the proper time formula (125) for G F which leads to: 



}jd 3 p e 




+C ° dT e e iT ^- m * +l ^T\S I y Pe~ Y Pe\- (129) 



Now, as this includes the momenta conservation constraints, this amplitude can be re-casted in 
terms of the dual (face) variables and thus easily written in terms of spin foam amplitudes. 

In the Euclidean context considered up to now, we ave been working with the Hadamard prop- 
agators. In order to properly recover quantum field theory amplitudes, one needs to construct 
(particle) observables which would reproduce the Wightman propagators and the Feynman propa- 
gators. These notions being closely related to causality are more likely to be properly implemented 
in a Lorentzian model only. The Lorentzian analog of the Ponzano-Regge model has been con- 
sidered in [26] and is known to be based on the group SU(1, 1). Therefore, we will study in the 
following how to write a (well-defined) Lorentzian spin foam model based on SU(1, 1), its abelian 
limit as a spin foam model over the abelian group M3 (Minkowski space) and causality oriented 
particle insertions. Discussion of implementation of causality in spin foam models can be found in 
[27]. 

An interesting feature of the Lorentz group SU(1, 1) is that there already exists a split between 
some so-called future oriented timelike representations and past oriented timelike representations. 
Moreover, the corresponding characters is then ex.p(=fidj9) instead of sm(dj6): their abelian limit 
will directly be the Wightman functions instead of the Hadamard functions. As we will see below, 
there is also a very natural implementation of the Feynman propagator. 

7.2 The abelian limit of the SU(1, 1) characters 

SU(1, 1) has two Cartan subgroups. The first one corresponds the (space) rotation and its elements 
are of the form (0 < 9 < 2tt): 

cos 9 sin t 
- sin 9 cos t 



hg 
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while the second one is made of the boosts (t € R) 



-t 



Then the (unitary) representations of SU(1, 1) are also of two kinds. We first have the continuous 
series of representations, labelled by s G R. They correspond to the (co-adjoint) orbits defined by 
the one-sheet hyperboloid in Kirillov's theory, The characters are: 

Xs {h e ) = Xs(±<h) 



sinht| 



Then we have two discrete series of representations -positive and negative series- whose correspond- 
ing orbits are the upper and lower hyperboloid. We label them by an integer n > 1. The characters 
are: 

Te ±i(n-1)6 e -(n-l)|t| e -(n-l)l*l 

XM= 2isin9 X « W = ^ X ^- at) = { - ir ^nhI\- 

Now let us consider the abelian/classical/flat limit of the characters for the (positive) discrete 
representations. We take the angles = itik (or t = run) going to zero while we take the represen- 
tation labels (n — 1) = 1/k to oo by taking the limit k — > 0. In this abelian limit, group elements of 
SU(1, 1) are mapped into Lie algebra vectors of su(l, 1) ~ M^y. rotations correspond to time-like 
vectors and boosts to space-like vectors. The limit of the character reads: 

_ e il\x\ 

Xn( h e) Xl(x) = for x timelike, (130) 

and 

e -l\*\ 

Xni. a t) Xi( x ) = ~^~T for % spacelike. (131) 

Then from the point of view of the abelian limit of the Lorentzian model, it is natural to consider 
Feynman evaluation with propagators given by G m (x) = Xm{x). Using the Kirillov formula: 

Vz G M 3 , G m {x) = m f d 2 ne im3 ^ = — [ d 3 p5(\p\ - m)e^' ', (132) 
Jn+ m Jc+ 

where is the future light cone and C+ the future cone. Comparing this expression with the 
definition (119), it is obvious that the propagator defined as the abelian limit of the character Xn ^ s 
the positive Wightman (Green) function. Similarly, if we would assign a negative representation n~ 
to an edge of the graph, we will put assign the negative Wightman function to that link. Finally, 
we could assign the representation n + © n~ to the edge, the corresponding character would be the 
sum xt, + Xn which would lead to the usual Hadamard (Green) function. 



7.3 Constructing a Lorentzian spin foam model 

Let us start by constructing an abelian Lorentzian spin foam model, with observables for particle in- 
sertion which exactly reproduce the Feynman amplitudes. Then we would like to write a Lorentzian 
version of the Ponzano-Regge model, based on SU(1, 1) whose abelian limit (corresponding to the 
Newton constant going to 0) reduces to that model. 
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Formally, we can define the abelian Lorentzian spin foam exactly the same way than the Eu- 
clidean one: 



Z = T\[ d 3 u s TT <$(«/), (133) 

where we now label the (dual) edges by vectors in the 3-dimensional Minkowski space M3. To 
insert a particle on an edge e, we identify the momentum p e with the geometric variable uj living 
on the dual plaquette. Then: 

• to insert the Hadamard function: one replaces the constraint 5(u e ) by 

5(\u e \ 2 -m 2 ). 

• to insert a Wightman function, one splits the <5(|tt e | 2 — m 2 ) in two depending on the time 
orientation and uses 

e(eu° e ) 5{\u f \ 2 -m 2 ). 



to insert the Feynman propagator, one splits <5(|il e | 2 — in 2 ) in two in the proper time formu- 
lation and uses 

p+00 

/ dT e i T e{\ue\ 2 -ml+ia) 

Jo 

where a is just a regularization. 



Now we would like to introduce the corresponding propagators in the non-abelian case, i.e. 
the observables describing particles in the Lorentzian spin foam model. Technically, our aim is to 
identify the quantity " A(9)8g(g)" which we should substitute to 8(g) in the spin foam model. Our 
main criteria is to check that the non-abelian propagators have the right expected behavior in the 
no-gravity limit n — > 0. 

Let us start by constructing the Hadamard propagator. The non-abelian counterpart of <5(|u e | 2 — 
m 2 ) is to fix the group element g e in the same conjugacy class than the Cartan element hg. Thus 
we introduce the distribution 69(g) defined by: 

/ dg5 e (g)f(g)= [ duf^hgu' 1 ), V/. (134) 

JSU(1,1) ^SU(l,l)/U(l) 

Let us insist at this point that h$, 6 G]0,7r] is not conjugated to h-g. Then as SU(1,1)/U(1) ~ 
TL + U TL- is the (disjoint) union of the upper and lower mass-shell hyperboloid, we can actually 
write: 

/ dg5 e>0 (g)f(g) = duf(uh e u- v ), / dg 5_ e (g)f(g) = duf(uh e u~ l ). 
JSU(1,1) JH+ JSU(1,1) J -ti- 

lt is easy to check that 

5 9 (g) =Y,Xn e (h9)x € n(9), (135) 

since the character of the continuous representations vanish on the U(l) Cartan elements hg. The 
Hadamard propagator is then defined as: 

s Hadamard {g) = 5(){g) + 5 _ e{g) . (1 36 ) 
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In the abelian limit, 9 G [0, tt] will be put it k units, 9 = ran, with m the (renormalized) mass. 
Notice that xt (~ ®) = Xn W> so that we can write: 

6 { e H) (9) = Y,XnWXn(9)- (137) 

n,e,e' 

Then the Wightman propagators only impose p° > or p° < 0. This is exactly whether u G 
or ii E K_, i.e. distinguishing 9 > from 6* < 0. Thus we define the non-abelian Wightman 
propagators for 9 G]0,7r]: 

4 W+) (5) = 5,( 5 ), 5f-\g) ee 5_ e {g). (138) 

The hard part is to define the non-abelian equivalent of the Feynman propagator. Looking 
at the proper time expression (125) of the propagators, we can interpret (137) as the equivalent 
expression of the Hadamard Green function in the non-abelian case. Then the Feynman propagator 
should be a particular splitting of the sum over e's. To decide, let us look at the abelian limit of 
Yl n Xn(@)Xn(9)- Using the explicit formula for the characters, it is easy to get (\p\ G 

£ mx £ n (0)xi(9) kZq J q +0 ° dl ^^e^)^), (139) 

when g = cos <p + is g(po) sin (pu. a, (f) G [0, 7r], u normalized positive timelike vector, is conjugated to 
the compact Cartan group. Here we have defined up = sg(po) sincfru for k — > 0. If g is conjugated 
to the non-compact Cartan group, it corresponds to a spacelike momentum and the limit is (still 
\p\ G R+): 



A(9) X l(9) X i(g) ~ / dl-^-e- l ^ m \ 



Inserting a regulator a > in the expression (139), we get: 

Y j e- (n - 1)a A(9)xmX € n(g) ~ / dl e i«e's 9 (ro)\p\+em + i a)j 
„ Jo 4 \P\ 



ee' s i 



»c-»o 4 |p| e'sjpj + em + ia ' 



(140) 



where we note s = sg(po). 

A priori, one would say that combining these expression, one can obtain all the different prop- 
agators of quantum field theory with different choices and signs of the poles. It is actually true 
except for the case of the Feynman propagator whose case is slightly more subtle. 

Indeed, we can combine 

1 1 1 

+ -r, = Hp\ 



s\p\ — m + ia s\p\ + m + ia \p\ 2 — m? + isa ' 

which would give the Feynman propagator if sa > 0. Therefore the sign of the regulator would 
need to depend on the sign s of the group element g. The easiest way out is to introduce a group 
element \g\ defined by: 

g = uh^u" 1 = cos <pld + i smcfiu.a =4> \g\ = uh^u^ 1 = cos <pld + i\ sin <p\u.a, 
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so that always s(\g\) = +1. Then it is straightforward to check that a good choice of Feynman 
propagator is: 

5f\g) ee 2iJ2e-^ a A(e) X -(e) X n(\9\) = 2i £ X n(M)(A(0 - ia) X -(0 - ta)). (141) 

n n 

since its classical limit is the expected l/(\p\ 2 — m 2 + iot). 

Conclusion 

In this paper we have shown how the quantum gravity amplitudes for particles coupled to three 
dimensional gravity proposed in [3] are related to the usual amplitude in the no-gravity limit. 
More surprisingly, it turns out that these amplitudes can be understood as coming from a non- 
commutative braided quantum field theory. 

We have seen however that the natural quantum gravity amplitudes amount to insert Hadamard 
propagators, then we have discussed how the insertion of a Feynman propagator can be naturally 
implemented in the spin foam model at least in the Lorentzian case. What would need to be done 
in this context is to show that this procedure is really well-defined in the spin foam context, namely 
that our proposal for introducing the Feynman propagator in the Lorentzian Ponzano-Regge model 
respect the topological invariance of the theory and leads to triangulation independent amplitudes. 
Also it would be interesting to understand more generally the definition of the Feynman propagator 
in the case of non-zero cosmological constant involving quantum group. 

It is not totally surprising that the effective field theory describing the dynamics of particles 
coupled to three space-time dimensions is given in term of a non-commutative braided quantum 
field theory. Indeed it is known that particles in 3d gravity respect the principle of Doubly special 
relativity and that the corresponding algebra of symmetry is a K-deformation of the Poincare group. 
Quantum field theories based on K-deformed Poincare group have been already considered in the 
literature [28]. Usually there is an ambiguity in the definition of the spacetime coordinates which 
leads to different proposals, none of which are physically motivated. Our analysis gives a unique 
prescription and solves this ambiguity. It seems that the explicit form of the action and star product 
which we found has never been proposed. An interesting perspective appear if one remark that 
the structure of the effective field theory which arise here possess some strong similarity with the 
Group field theories that appears in the spin foam quantization of gravity [30], it is tempting to 
conjecture that this similarity is not purely accidental and that the structure which is unravelled 
here should generally appear in the spin foam context. 

This questions becomes more interesting if we think of 2+1 gravity as a toy model for 3+1 gravity 
and hope to learn something about the effective dynamics of matter coupled to quantum gravity. 
This hope is substantiated by the recent work [29] which shows a reformulation of 3+1 gravity in 
terms of a perturbed BF theory where matter coupling to gravity is very similar to the 2+1 gravity 
case. Moreover, it seems that the star product that we have introduced should be generalized to 
a 3+1 quantum space-time. A subtle point is that star products in the field theory action would 
then also explicitly involve polynomial terms in the coordinates (or equivalently derivations in 
the momenta), which are usually ignored in standard analysis of possible non-commutative field 
theories as effective quantum gravity theories. This specific structure allows the field theory to be 
actually invariant under the deformation of the Poincare group. 
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Figure 4: Tetrahedral net as the skeleton of a 2-sphere. 
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A Examples of the duality transform of Feynman diagrams 

A.l Example of T ^ S 2 

Let us consider the tetranet embedded in S 2 . We have six edges thus six (momentum) variables 
o, b, c, d,e, f, and four vertices resulting in four constraints: 

' a+b+c=0 

a — d + e = 

c+ d+f=0 
, -b + e + f = 

The rank of this system is 3, so that we need at least 6 — 3 variables to parametrize its solutions. 
A convenient parametrization is given by assigning a new variable to each face, thus four variables 

a = a — 5, b = (3 — a, c = S — (3, d = 7 — 5, e = 7 — a, f = (3 — 7. (143) 

It is obvious that we can translate all face variables by a constant vector k and it wouldn't change 
anything. So, at the end of the day, we truly have 3 parameters to describe the space of solutions, 
which is consistent. 

More formally, one can introduce conjugate variables to the face variables a, (3, .... Let us note 
them v a , .... Then C = v a + vp + + v$ generates simultaneous translations on a, so that we 
can see the residual invariance as imposed by the constraint (7 = 0. This can be understood as the 
closure constraints for the tetrahedron assuming that the vectors v a , ... are the normal vectors to 
the faces of the tetrahedron. 

A. 2 Example of T ^ T 2 

One can also embed the tetranet in the torus T 2 . As the graph is the same as in the previous case, 
we have the same variables with the same constraints. The difference is that we now have only two 
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(142) 




Figure 5: Tetrahedral net as the skeleton of a 2-torus. 



faces a, (3. This is not enough variables to parametrize the space of solutions. We need two extra 
vectors, which will be associated to the two (non-contractible) cycles of the torus: 

a = d 

b = (a - (3) - d 

c = {(3 — a) 

d = (a - (3) - C 2 

e = (a-(3)-C 1 -C 2 

{ f = c 2 



< 



The chosen cycles, as sequences of edges, are a,—b,—e and —d,—e,f. Let us also notice that 
once again we can translate a, (3 simultaneously by any vector without changing anything. This 
corresponds to the closure of the faces. 



A.3 Example of Y ^ S 2 ^ S 3 

Let us consider the 3d manifold 5 3 and triangulate it using two tetrahedra. The skeleton of the 
resulting triangulation A is simply a tetranet, which we take as on fig. 4. Let us then consider the 
Feynman diagram T made of the edges a, b, c, d, e. The graph has four vertices and the corresponding 
equations are: 

b = e, c + d = 0, a + 6 + c = 0, a + e = d. (144) 

To parametrize the space of solutions, we introduce embed T in the tetranet considered as a 
triangulation of S 2 . The tetranet has four faces for each of which we introduce a new variable: 
u\ for abe, u 2 for acd, 113 for bcf and U4 for def. Then we can express the original momentum 
variables a, .., e in terms of these face variables: 

a = u\ — u 2 , b = U3 — ui, c = u 2 — us, d = 114 — u 2 , e = U4 — u±. 

For this to be a solution of the constraints (144), it is necessary and sufficient to impose / = 
us — «4 = 0. And we finally recover the formula (10) for the Feynman amplitude. 



B Feynman graph with Spinning particles 
B.l Overview of the spinning particle 

In [3] Ponzano-Regge amplitudes for spinning and spinless particles have been proposed. In the 
main text we have shown that the spinless Ponzano-Regge amplitudes reproduces in the no gravity 
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limit the usual Hadamard-Feynman amplitudes. In this appendix we show that the same is true 
for the spinning amplitudes. Let us present a quick review of the properties and propagator for 
a massive spinning particle/field in three dimensions. Let us start by the Euclidean case, keeping 
that all the considerations could be straightforwardly translated to the Minkowskian case. 

Let us look at a one-particle state ip PjCT where p denotes the (3-)momentum and a the other 
degrees of freedom: 

P^p,a = P^p,a- 

Considering the action of Lorentz transformations U(A) on the translation generators, U(A) maps 
a state of momentum p onto a state of momentum Ap: 

U(A)ip Pj<T = ^C CT >(A,p)V>A P , CT '- 

a' 

Let now us choose a reference momentum k^ and define states ip Pja from the states ipk,a : 

where N(p) is a normalization factor and s(p) maps k to p. For a massive particle, k is usually 
taken as p^ ^ = (1,0,0). The little group of Lorentz transformations W leaving k invariant is the 
rotation group U(l). So s(p) is a section 16 from SU(2)/U(1) to SU(2) such that: 

s(p).(l,0,0) = -p= (Jl-^,-p). 

m \ V m 2 m I 

Then the action of a Lorentz transformation reads: 

[/(A)VWx = N{p)U(s(Ap))U{W{A,p))^ a , 

with W(A,p) = s(Ap)~ 1 As(p) living in the little group U(l). Note that W(R,p^) = R for rotations 
and W{B,pW) = Id for pure boosts. From this, it is straightforward to check that postulating the 
action of the little group: 

a' 

we get the action of an arbitrary Lorentz transformations: 

U{A) ^ = (iv^)) T.B'Aw&pMaps- 

So that instead of giving all the coefficients C CT / )(T (A), it is enough to define the coefficient D c i^{W): 
from a representation of the little group U(l), we induce a representation of the Lorentz group 
SU(2). The usual normalization is N(p) = y^o/Po so that the scalar product is normalized 
(i^p,cr,^p,cr) = S aa rS^ 2 \p — p') . The spin s is the choice of a irreducible representation of the little 
group U(l), which defines the coefficient D a i <a , and a is the angular momentum. 

16 We choose 

s(m cosh #,msinh# cos <?i, msinh#sin0) — e <t " T3 e 6a2 . (145) 
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All this can be directly translated to a field theory. Let us assume that the field ip has compo- 
nents ipi labelled by I and is a representation of the Poincare group. More precisely, we split the 
field into creation and annihilation parts tpi = ^ + + V'; - : 

tt = E f*?«M <*<* *r = E / *p«m <*.'<* 

a a 

Then the action of a Poincare transformation reads: 

U{K,a)^fU{K,a)- 1 = £ D^A'^iAx + a), 

V 

where we have chosen a representation I os the Lorentz group SU(2). Interactions for such a field 
can be written as: 

V(x) = g h ...i N r 1 ...v M ^ i {x) • • • Y^(x)Y>+(x) . . . ^ N (x), 

where g^... is an intertwiner for the Lorentz group (in our case where the Lorentz group is SU(2), 
it is a Clebsh-Gordan coefficient), so that the interaction is a scalar for the Lorentz group. Then 
the coefficients u\ (and v\) defining the field satisfy the following constraints: 

u l (x,p,a) = (2ir)- 3 / 2 e i v x u l (p,cj), 



^u v ((Ap),a')D s a , a (W(A,p)) = ^EA'i(A)«i(p». 

In order to understand the structure of this constraint, let us split it in the two cases when A is 
a pure boost or a rotation. A boost, with W(A,p) = 1, allows to map the reference momentum 
( 2 )p = to an arbitrary momentum: 



ui'{p,cr) = J —^2Di>i(s(p))ui(0,a), 
V P l 

so that we go from ui(0,a) to ui(p,a) simply by a boost. Then for a rotation A = R, with 
W(R,p) = R and ^p = (Rp) = 0, we have: 

]T u ,,(0, v')D s a , a {R) = £ L^(i?H(0, a), 

cr' I 

which states that the coefficient ui(0, a) can be understood as an intertwiner between the represen- 
tation of U(l) defined by the spin s and the representation of the little group U(l) induced by the 
representation I of the Lorentz group SU(2). So that the field definitively lives in the representation 
of spin s of the little group. In our Euclidean 3d context, I is a representation of SU(2) and s, the 
spin, a representation of U(l) such that s < I (i.e s <^-> I). There is only one value of a which is s 
and D S (R(9)) = exp(is9) for a rotation of angle 9. Then the unique solution to the constraints is 
«,(0) = (l\s) = Su. 

Finally, the causal propagator is defined as the commutator of the field ip with itself. Assuming 
that vi = uj, we have: 

W(x)Mv)] = (2^)3 / d 2 pA ul (p)e^-y\ with A w (p) = £ a)u\,{p, a). 
To sum up, the propagator A^,(p) of a field of spin s must satisfy the following properties: 
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• A(ApA _1 ) = AA(p)A- 1 . 

• A(p(°)) is the projector \s)(s\ to the right spin. 

• A(p) is a projector: A(p) = A(p)A(p). 

I represents the total angular momentum while s, the spin, is only the intrinsic rotation of the 
particle/field: \s) is a vector of the / representation i.e s < I. From the field point of view, I = s 
represents the fundamental field 4>^ of spin s, while the representation I = s + k corresponds to 
the k-th derivative of the field d k cj) <yS \ 

The unique solution to these constraints is: 

A(p)=D I (s(p))\s)(s\D I (s(p))- 1 . (146) 

This is easy to generalize to different Fs on the left and right hand side. Now 

A(p) = D I (s(p))\s)(s\D 1 ' (s(p))- 1 

corresponds to the correlation (d k <j)^ d k '(j)^) with I = s + k, I' = s + k' . 

For I = s = 0, the scalar field, the propagator is obviously A = 1. For the spin-1/2 particle, 
I = s = 1/2, the definition of the section s(p) reads in spinor notations 17 : 

where the dj's are the Pauli matrices normalized to o; L Oi = Id. So the propagator for s = +1/2 is: 
A+i /2 (p) = S*(*(P))(J l^DHslp))- 1 (147) 
= DH*(p))\(ld+( I ^ 1 ))l>'(*(p))- 1 = ^(m+^^) > 
and similarly for the negative helicity s = —1/2: 



(148) 



A+i/ 2 (p) = ^ (m -!?■<?). (149) 

And we recognize the usual expressions for a spinor field for the two helicities. 
For a / = 1 field, the section is defined by: 

D^isip)). f )=-[ Z ) => D(s{p)) = R{p)B{\ 2 p\)R{Pr\ 







in 




7 The spinor representation of a 3- vector (110,11!,%) is: 



M() u + = u x + iu y 

U- = u x — iuy —uo 
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If we look at the spin s = field, the projector is: 

|/=1,0)(0| = 

so that the propagator is 




7=1, s =o 



(p) = 



which naturally corresponds to the correlation (d^(pd u (p) of a scalar field <p. To recover the usual 
propagator for a spin-1 field, s = 1, we in fact sum over the states |s = —1) and |s = +1). The 
projector then reads: 



-X-l + I+X+I 




so that the propagator is simply: 



A 



1=6 



"\v) = V - 



(150) 



B.2 Inserting spinning particles in the spin foam 

Now, let us look at the Feynman amplitude of a graph with massive and spinning particles. The 
part concerning the mass will not change: we still impose the constraint 5(\p\ — m) on each edge 
and conservation of the momentum at each vertex. Moreover we now need to take into account 
the further degrees of freedom corresponding to the spin (or intrinsic rotation) of the particles. 
On each edge e, we have a spin s e and we put the corresponding propagator A /=Se (p e ). At the 
vertices, we need to take into account the interaction between spinning particles: we include SU(2) 
intertwiners intertwining between the representations / = s e of the edges meeting at the vertices. 
On the whole, the Feynman evaluation reads: 

* = / n ^ n '"";!,; ""V r w * n 5 f e a - e A °t, (m 

e e e v \e\v=t(e) e\v=s(e) j 

Let us point out that we need not in principle take I e = s e , then the propagators would live in 
the representations I e and the intertwiners C intertwine between the I e 's: the spins s e would only 
enter as projectors in the definition of the propagators A^/(p) = Di s (s(p))D s ii(s(p))~ 1 . 

It is straightforward to write this amplitude in its dual form solving the constraints imposed 
by momentum conservation at the vertices. Indeed, adding spins doesn't modify anything, so that 
the amplitude corresponding to the graph embedded in the 3 triangulation reads: 

Ir=f n^e-II*K)II 6{W t^ me) A t^(Pe) X Ct- (152) 
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So given the graph T, we identify the plaquettes dual to the edge of T. These form a tube around 
the graph. We evaluate the holonomy p e around each plaquette, then joining them all into a 
spin network -a graph drawn on the surface of the tube- using the intertwiners C\ ". Finally the 
Feynman evaluation id the product of the mass shell conditions with the spin network evaluation. 

Let us now transpose this formula to the non-abelian case in order to encode spinning particles in 
the Ponzano-Regge model, such that its abelian limit will reduce to the above Feynman amplitude. 

The mass shell condition translates into Sg e {g e ), which imposes that g e = uhgu^ 1 for some 
u £ SU(2)/U(1). u represents the momentum of the particle so that the spinning propagator 
translates to D[^ =s (u)D I s jT s (u~' 1 ). Finally, we contract these propagators using the intertwiners C's 
which defines the physical content of the quantum field theory we would like to encode on the spin 
foam. The final amplitude is: 

Zr, { e eM = f II II pjfcn L*mm ^^""X^ H<T (O * ( 

J g e ^p eer <j \ ej JbV(2)/U(l) 

(153) 
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